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EXECUTIVE SUMMARY 

In current missile designs, it is common practice to form the 
main propulsive system by clustering po,gether several rocket engines. ' 
The clustered engine concept is advantageous to many missile design 
dis ciplines. in that it tends to. optimhze the vehicle. However, this design 
creates a very complex thermal fluid problem commonly referred to 
as the ”3ase Heating- Problem". Heating of the base region results 
from radiation which is emitted by the high-temperature rocket exhausts 
and-'convection from, gases circulating in the base region. As a result, 
a Severe thermal environment usually exists in the base region. Vital 
engine components and sensitive base structure must be protected from 
this environment. This is accomiplished by determining the thermal 
response of the components to the' base environment and designing an 
adequate protection system. Since this must be accomplished prior to 
any flight tests, the ability to predict the environment is a necessity. 

As a result of experience gained in the design of previous rocket 
vehicles, several methods of predicting the base region thermal environ- 
m.ent have evolved. , These can' be generally categorized as: 

G -Scale model testing 

o Extrapolation of previous and related flight test 
results. 

o Semiem.pirical analytical techniques. 

A study of both the theoretical and experimental aspects of the problem 
will show, however, that all techniques are still in the development stage 

A. NASA-funded review of the base heating experienced by the 
Saturn Apollo vehicles was initiated in July 1971. The objectives of 
this - study were the following: 


4 ; 


o Collect, review, correlate, and compare all relevant 
Saturn flight test, scale model test data, and analyti- 
cally produced base environmental paranieters 

o Identify, ^ delineate, and e:mlain the deficiencies and • 
differences in the experimental and analytical 
methodologies 

o Develop improved analytical and experimental pro- 
cedures for predicting and evaluating the base . 
environment. 


Results, conclusions, and r ecomnicndations derived from the study are 
contained in three volumes. The accomplishments and conclusions 
resulting from this study are summarized in the following paragraphs. 

Pertinent Saturn vehicle flight and scale model test data were 


collected, correlated; compared, and are presented in the reports. 

The review of the scale model test results indicated the 
fcllov/ing: 

o Scale model testing of the type used during the Saturn • 
prograna does not predict the full scale vehicle radiant 
environment (Note; The Saturn model test program 
\vas> designed only to provide information- about the 
convective environment. ) 

o Scale model testing tends to over-predict the con- 
vective environment • ‘ . 

c Test results need to be scaled or otherwise altered 
in developing full-scale base thermal design criteria. 

At present, the scaling procedures produce uncertain 
predictions of full scale results. 

o Scale model test data have been used with better 
results in establishing the design criteria for upper 
stages than for booster stages. 


o Scale naodel tests have been found useful in establish- 
ing the relative effects of some of the off-nominal 
vehicle operating conditions (Note: Exceptions to this 
can be fo'ond. ). 


I 






and o 

base 

ance 


It was concluded th.at although scale model tests are expcnsiv^c 
ften leave much to be desired in establishing the iviil-scale vehicle 
thernaal design criteria, if properly used, they can provide guid- 
and a degree of confidence in the tecimique which is being used. 


Work associated with the analytical methodologies produced 
the following: 

p A procedure for predicting the convective cooling 

during the aspiration phase of flight was developed. . 
Comparisons with the flight test results are pre- 
sented in the report. 

o The computer technique .for predicting the base 
environment during flow reversal was improved, 
expanded, and computed results compared to. 
measured flight and scale model test data. ‘ The • 
comparisons which are presented in the report gen- 
erally indicate that .the progra'm tends to conserva- 
tively over-predict the base convective flux, over- 
predict the gas temperature, and under-predict the 
pressure. Coupled v/ith the comparisons contained, 
in the -report, the flow reversal program should'pro-' 
vide an alternate method for developing full-scale, 
vehicle base design criteria. 

o The use of NASA- developed gaseous radiation com- 
puter programs appears to offer an appropriate- 
method for predicting plume radiation^ However, 
this method is somewhat complicated and time-, 
consuming. 

o A simplified radiation prediction technique used . 

■ successfully to predict the Saturn and Titan base 
radiation is easier to use, but requires a knov/1- 
edg^e of the plwme ground-level e-missive power. A 
cornputer program was developed to compute the 
radiant fliu-x in the base region using measured plume 
emissive power as input. Partially occluded plumes 
and blockage by nonparticipating surfaces can be 

. incruded. The change in radiation v/ith altitude is 
obtained by extrapolating results measured during 

■ flight tests of sinailar vehicles, which is one of the 
li'rnitations with this technique. 


o 


Stucaes of cr.yogcnically cooled, eagir.es aidicatcd ■ 
that condensation of the exhaust gases along the 
engine wail can occur when the wall tcraperature . 
is beiov/ the saturation temperature of the exhaust 
gases. The effects of condensation on tiie base 
environment could be significant. 


Review of the flight test results produced fee following: 


o Useful qualitative and quantitative trends and 

characteristics can be obtained from the measured 
flight test data. 

o All instruments used to measure the Saturn inflight 
•base -heat transfer parameters had the potential for 
erroneous indications. The quantitative accuracy of 
some of the instruments could not be determined with 
the information available. 


o Techniques were developed for correcting some of the 
measured results. For future flight tests it is recom- 
mended that the instrunaents be calibrated in an environ- 
ment similar to that anticipated during the flight, i. e. , 
a combined radiant and convective environment. The 
- necessary flight parameters can then be evaluated, 
approximately, from the calibration data. 

The procedures and methodologies for predicting the base 
environment and developing thermal design criteria are basically 
developed; however, many of the complexities of the problem will 
never be significantly simplified. Many improvements can be made 
and further research is recommended. Areas and specific problems 
which should be studied are listed at the end of each section in Volume I 
of this report. 
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M. INTRODUCTIOiM AND SUfvlMARY '.. 


1. 1 


INTRODUCTION 


In current mis sile designs, if is common practice to form the 
main propulsive system by clustering together several rocket engines. 
Figure 1-i shows two typical ar r angements which have been employed. 
The clustered engine concept is advantageous to many missile design'- 
disciplines,, in that it tends to optimize the vehicle by resulting in a 
relatively compact, minimum weight, and dynamically stable system.- 
'However, this design creates a very complex thermal-fluid problem 
commonly referred to- as the ’‘Base Heating Problem”. The heating 
results from the radiation emitted by high-temperature rocket exhausts 
and convectively from gases circulating in the base region. This phe- 
nomenon is shown schematically in Figure 1-2. In summary,' the heat- 
ing in the base region varies significantly with altitude- and can be 
characterized as follows: at low altitudes, the supersonic exhausts, in 
conjunction with the free stream flow, pulls' the ambient air into the base 
region which provides cooling to offset the usually high thermal radia- 
tion. Ar high alnt'udes, where the ambient. pressure is low, the exhaust 
jets expand and interaction between adjacent engine plumes occurs. 

This causes some of the high-ternperature exhaust gases to reverse 
back into the base region, resulting in convective heating. Generally, 
heating rates and gas temperatures in the base region are relatively 
high, and fluxes of 10 to 30 Btu/ft^-sec and temperatures of 2, 000®F 
are not uncommon. 


As a result of the severe thermal environment, . vital engine com- 
ponents and sensitive base structure must be protected. This is accom- 
plished by determining the thermal response of the components to the 
base thermal environment and designing an adequate protection system. 


1-1 
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FIGURE 1-2. BASE FLOW RE6IM 



bince *chis must be accompli shec prior, to any flight tests, the ability t*o 
predict the enmrcnment is a necessity. Unfortunately, because of the 
very complex narure of the base heating problem, the pr cdictability of 
the en^h.rorimer:r is very limited, even though much time and liioncy 
have been devotee to researching the various aspects of the problem 
during the p>ast.lO years. 


As a. result of previous efforts, several methods of predicting 
the. base region thermal environment have evolved. These methods 
can be generally categorized as: 1) Scale model testing; 2) Extrapolation 

of premous and related flight test results; 3) Senai -empi rical analytical 
techniques. ■ 


However, in studying both the theoretical and experimental 
aspects of the problem, it becomes obvious that all techniques for pre- 
dicting the base heat transfer are still in the development stage. This 
is not surprising when the extremely complex nature of the problem is 
considered. The theories associated with such problems as fully 
separated and interacting three-dimensional compressible viscous 
fluid flows, multiphase turbulent and chemically reacting boundary 
layers and coupled radiation, and convection and conduction heat trans- 
fer (to mention just a few) are currently limited and complicated. How- 
ever, these problems are commonly encountered in the base heating 
problem. In the development of the analytical techniques, these com- . 
piex problems are solved using much simpler theories which ultimately 
introduce inaccuracies into the calculated results. The basic idea in 
the analytical model development is to approximate the insolvable pro- 
blems by simple theories and then to adjust the computed results until 
agreement w^ith experimental results is obtained. Thus, the. final ana- 
lytical model would be a semi -empirical solution to the problem. 
Trouble is encountered using this approach when considering the 
experiir.entai data Vv^hich is to be used as the basis for modifying the. 



ar.ixiytical model. If scale -model. te st .results are used, these results 
often do not agree with the flight test results because of the inability 
to simulate full-scale vehicle conditions. The flight test results, 
on the other hand, are suspected of being in error at various times 
curing the flights. As a result, if the analytical .models are forced 
to a.gree with the available experimental data, a general solution of 
the base heating problem is not obtained. IVhat does emerge is a 
model which fails to properly predict the environment for vehicle con- 
figurations which are significantly dissimilar. 

Other factors also contribute to the general lack of base heat 
•transfer predictability. Such realistic factors as multiphase (i. e. , 
engine wall condensation) engine boundary layers, .surfaces with tem- 
peratures controlled by cryogenic liquids, and water absorption by 
porous heat shield materials appear to affect the base environment 
directly. These factors are not the same on all vehicles, and the 
predictability of the effect is often speculative. 

1. 2 SUMMARY’ ■ ■ ■ ' 

In this, study all available analytical techniques, Saturn flight 

i 

test results and scale -model test data were to be reviewed. The 
purpose of this review was to determine how v/ell these results a.gree 

to provide guidance in establishing acceptable techniques for 
ctinr/ the base thermal environment and base tnermal protection 


r— i 




o Task 1 - Fvevi.ev/ of Analytical Methods and Svaluation of 
IrnAt^ortant parameters - All available analytical methods 
(for both solid and liquid propellant systems) for caL ula- 
tion of base heating data were reviewed and evaluated. 

The validity of assumptions made in the analytical methods 
was investigated. The most influential parameters were 
determined. Condensation, turbine exhaust gas injection, 
shroud, chemical kinetics, and other effects that may 

influence base heating were considered. 

» 

o Task 2 -^ Evaluation of Commuter Programs - Available 
commuter nrogram.s and mianuals for calculation of base 
'heating were reviewed. Representative analytical 
miethcds and Teledyne Brown Engineering's base heating 
m.anual were used to generate analytical' data for Saturn 
vehicles which were compared with flight and experi- 
mental data in Task 6. 

o Task 3 - Collection and Review of Experimental and Flight 
- Data - All avaiiabie Saturn I, IB, and V flight data, static 
test data, and scale miodel test data, including both hot and 
• cold flow tests, were collected and reviewed. These data*- 
were separated and categorized for the required correla- 
tions, scaling, and corrections w^hich were performed in 
Tasks 4 and 5. 
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Task 4 - Correlation of Flight and Experimental Data - 
Measured data associated with each stage (Saturn I, IB, 
and V) that best characterizes the base heating phenomenon 
were selected for evaluation. Off-nominal conditions were 
identified and investigated. Data correlation, correction, 
and comparison techniques, were studied. - Errors associ- 
ated with mieasurem.ents and instrumientation were investi- 
gated. The rniost suitable technique was used to correlate 
and compare the measured data. Factors which cause 
correlation deficiencies and which create off-nominal 
conditions were identified and investigated. 

Task 5 - Scaling of Model-Test Data - All existing scaling 
miethods \vere reviev/ed and evaluated. The deficiencies 
of existing scaling m.ethods v/ere determined and improve- 
mients made* The miodel test data were then scaled (by 
the best available scaling method) to full-scale flight 
dara. 
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Conclusions 


Recom- 



Task 6 - CouToarisons of Data, 
nocndations - The analytical data from' 


Task 2, the corre- 


lated flight and experimental data from Task 4, and. the 
scaled naodel test data from Task 5 v/ere compared and 
studi.ed. Base heating rates, recovery temperatures, base 
pressures, and other variables or phenomenon, e. g. , 
condensation effects, found to be important in the study 
were included in the comparisons* The objective of these 
comparisons is to explain the large differences between 
predicted and flight base heating data. 


o Task 7 - Preparation of Final Re~port - A final technical 
report that describes in detail the results, conclusions, 
and r ecomrnendations of. the base heating studies has been 
prepared. . • 


In Section 2, the basic theoretical aspects of the rocket base 
laeating problem are discussed. The fundamental equations, which show 
in simple termis what must be known in order to arrive; at an analytical 
solution of the base heating problem^, are discussed. 


In Section 3, the problem.s associated with the Instrumentation 
used to mieasure the inflight base heat transfer are discussed. The 
design of accurate instruments for measuring the base heat transfer is 
a difficult problem which has not been solved. -Many of the instruments 
used during the Saturn flight tests had the potential for error and appear 
to have produced mieasured results which cannot be completely corre- 
lated and evaluated. 


In Section 4, some typical Saturn flight test results ,are com- . 
pared, analyzed, and discussed. 

In Section 5, the scale miodel test results are compared to the 
measured flight data and discussed. 

Base heat transfer diagnostic and prediction techniques are dis- 
cussed in Section 6. Base heat transfer environmental parameters pre- 
dicted by somie of the analytical techniques are compared to the scale 
model and flight test results. 


During the course of the study, several quarterly progress 
reports were, published and are listed as References 1-1 through 1-6- 
Additional information concerning the Saturn base heat transfer can 
also be obtained from Reference 1-7. 
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2. DISCUSSiON OF THE BASE HEATING PROBLEM 


BASIC ?IEAT TRANSFER CONSIDERATIONS 


The heat transfer’ to the base of a clustered-engine rocket vehicle 
is a result of incident thermal radiation and convection combined with the 
conduction, heat storage, reradiation, and possibly phase changes asso- 
ciated with the.base region materials. Conside-r the differential material 
volume arbitrarily located in abase region as sho\vn in Figure 2-1. The 
differential area, dAg, is assumed to be exposed to the convective flux, 
and to the radiant flux, Assuming that the material is a porous 

mediuna, such as an insulation that could have absorbed moisture, the 
phase change which would take place as the material is heated at the 
saturation temperature represents an energy loss, q^. The reradiation, 
from the surface to the surroundings will also represent an energy 
loss. A simply energy balance' shows that the net heat crossing the unit 
differential area, dAs, is 


^in ” 9c ‘ ^9ir “ 91 9rr 


(2-1) 


net heat flux crossing dAs 


convective heat flux 


incident radiant flux 


- -reradiation flux 

- heat loss by phase change 

a - surface absorptivity. 

The temperature of the differential volume can be computed by solving 


Eouation 2-2, the Fourier conduction ecuation; 
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pcp = V(kVT) -f Qi 

*■ c • 


GX \; dX 


density 


temperarare 


thermal conductivity 


heat lost by phase change or gained because of 


'heat oreneration. 


Equation Z-1 is related to Equation 2-2 through the boundary 
conditions of the problem. For example, the conduction heat flvix is 
defined as 


^ - 1 II 

^ 3x 


(2-3) 


At the surface, dAs (see Figure 2-1), this is equal to the heat 
flux entering the Volume element and gives. the boundary condition; 


x=0 


x=0.- 


" °in 


(2-4) 


= qc' 4 G^q-ir - qi - qrr 


2-3 ■ 


1 ■ 
i ) 


^.'or^'laily Equation 2-E cannot be solved by closed-iorrn integra- 
tion for complicated heat transfer problems. As a result,- most problems 


are solved using a computer program such as the Chrysler Improved • 
Numerical Differe'-icing Analyzer (CINDA). This particular program'is 
very versatile and has the capabilities to include the effects of phase 
change heat losses, V/ith regard to Equations 2-2 and 2-4, the driving 
potential which causes the heat transfer is the convective and the radiant 
heat flux. As opposed- to this, in a transient problem, the conduction, 
reradiation, and the heat loss by phase change are controlled mostly by 
the internal properties of the material. With respect to base' heat trans- 
zer, it is usually the convective and radiant heating rates which are the 
most difficult to specify. • 
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BA.SIC PARACvlETERS IN CONVECTIVE HSA.T TRANSFER 


To show wr-'t must be known in order to specify the convective 
flux, consider the following. The convective flux is defined by Ne.wton's 
law of cooling as • , ' ■ ' 


q, = h (T^ - T^) 


(2-5) 


v/n ere 


hi 


w 


heat transfer coefficient 


jas recovery temperature 


surface temperature. 


Also, it is known that the convective coefficient can be derived using 
a modified form of the Reynolds analogy given as. 
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Nusselt n. umber = hx/k • • 

Reyr^olds number = purc/u 
Prandti number = ^c^/k 
experimentally determined coefficients 

characteristic dinnension 
thermal conductivity 
gas density 
velocity 

dynamic viscosity 
-* specific heat. 

Because convective heat transfer is a boundary layer problem, • 
property gradients exist between the outer fluid and that near the wall* 

As a result, the evaluation of the transport properties is dictated by 
correlating experimental data using Equation 2-6. The coefficients c, 
m, n are also determined as a result of this correlation. Thus, it can 
be said that' solutions to the convective heating are always semi-empirical. 

Substitution of the physical quantities into Equation2-6 gives 


u 


u 


V =: 


(c) • \ 


/ kV-^* p'^ 


a 


m-n 




(2-7) 
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Further simplification of Equation 2-7 is possible if the velocity, 
u, and the density, p, can be related to the local free&tream variables, 


. 2-5 
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1 1 


using ideal gas and isentronic relations. Then, 
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■ static pressure 

TO 
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total pressure 

2. g 

- 

static temperature 

V 

- 

ratio of specific heats 

R 

- 

gas constant 

M 

- 

local Mach number 

<y 

to 
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gravitational constant, 


- 1 


;RT. 
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Then, using the isentropic relation 

R /a° ) 


(2-10) 


where T is the total temperature. Equation 2-7 can be written as 
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( 2 - 11 ) 



which is 


constant 

terms 


/.transport 
\ properties 


/ characteristic \ / flow field 

\dimensicn / \ properties 

Equation 2-11 represents only one of the many relations which 
are possible. This equation relates the heat transfer coefficient to the 
most basic properties of the flow. . The computation of the heat transfer 
coefficient using Equation 2-11 requires the following knowledge: 

o The type, of fluid must be kno\%m in order to evaluate 
the transport properties Cp, k, |jl, y, and R 

o The flow field must either be known or be possible 
to compute in order to specify the thermodynamic 
properties Pg, Tq, and Tg 

o The coefficients c, m, n and characteristic dimen- 
sion, X, must be known or assumed. 

The latter implies that some geometrically similar problem must have 
been investigated and correlated using Equation 2-6. 
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With respect to trie base heating probienr, the following charac- 


u ! 


o Becaus e the base flow changes throughout the flight, 
all terms on the right side of Equation Z- 11 -are 
uransien'c 

o The gas in the base region during- the aspiration 
phase of the flight consists primarily of. the free- 
stream ambient air 


u 


I I ' 

U 1 


n i 


o During the flow reversal phase, . the base gases consist 
primarily of the reversed engine exhaust products 

o Regardless of how it is accomplished; a base gas . 

teinperatur e and a component surface temperature . ' > 

iTiUst be determined in order to solve Equation Z-1. 

Therefore, because the transport properties are primarily a 
function of temperature, it can be assumed that these properties will 
be known v/hen a solution is possible. Knowing this fact, the functional 
relations can be written 


,U ! 


I !•: 


qc(T) ■ =. £ [ h(T), Tr(T)'] 


(2-12) 


?he relations on the right side of Equation 2-13 are also' time-dependent. 


n i 


h(r) = f [ c, m, n, Pq, Pg, T^, x] 


(2-13) 


r**! i 


h(-) = f[c, m, n, Pg, Tg, u, x] 


LJ f 


D 


Tr(r) = f[T.(r)] 


(2-14) 


r~t 1 







Equations 2-12 through 2-14 are by no means a major simplifi- 
cation. However, these equations do show, in simrple terms, what must 
be known to derive an analytical solution* The analytical model must be 
devised so that the transient values of the pressure P^) and tem- 

perature (T^, T.^) can be evaluated from the known vehicle and thermo- 
dynaimic conditions. Also, the available experimental data must be 
correlated (using Eo^uation 2-6) to determine c, . m, n, and the chara;c- 


2.3 ASPIRATION PHASE CONVECTIVE COOLING 

After engine ignition and release’, the vehicle rises through the 
atmosphere with increasing velocity. During the, first phase of the 
ascent, the supersonic exhaust jets, in conjunction with the freestream, 
tend to create an ejector effect in the base region (see Figure 1-2), 

This action pulls the ambient air into the base region. The ambient 
air's having a relatively low total temperature during this phase of the 
flight provides a source of convective cooling to offset the usually high 
radiation. The aspiration phase usually persists only at low altitudes 
during the first half of the flight. Because of its positive effects, con- 
vective cooling has been largely ignored in base region thermal analyses. 
However, some attempts have been made to more effectively utilize this 
potential by placing scoops and fiov/ deflectors aro\ind the base of the 
vehicle (see Figure 1-1). The development of an engineering model of 
the convective cooling during this phase will be discussed in another 
section of this document. 


2.4 FLOW P.EVERSAL PHASE CONVECTIVE HEATING 

At high altitudes, where the ambient pressure is low, the exhaust 
jets expand until interactions between adjacent engine plumes occur. 

As a result of these interactions a very cornplex shock structure and flow 


< — ; 


iieicl is established (see Figure 2-2). Briefly, what happens is the 


ioiiov/in<^. 


n i' 


As the supersonic jet leaves the nozzle, turbulent mixing between 
the jet-free boundary and the surrounding gas takes place. As the jet 
approaches the impingement point, the mixing zone has spread to a rcla- , 
lively t'hick layer. Because of the mom-entum exchange associated with 
the mixing porcess, the velocity (and,, therefore, the streamwise momen- 
tum) varies significantly across the mixing layer. As the flow approaches 
the impingem.ent point, some of the low-energy gases in the mixing layer 
cannot penetrate the oblique shock wave and are therefore reversed. 

The criteria which establishes the portion of the mixing. .layer 
which is reversed is to determine the minimum energy streamline (i. e. , . 
the discriminating "D-'' stream.line in Figure 2-2) -which can penetrate 
the shock. , Tliis can be established by finding the streamline with a 
stagnation pressure equal to the static pressure behind the oblique 
shock. Likewise, the location of the undisturbed jet boundary (i. e. , 

”J" streamline in Figure 2-2) determines if any of the main jet. is 
reversed. If the " J" streamline of Figure 2-2 has a stagnation pressure 
greater than that of the discriminating ''D'' streamline, then the net 
flow is but of the base region (i. e. , no flow is reversed). Plowever, 
when the "J” streaimline stagnation pressure is less than that of the 
:iD" str-ea’Tiline, tlnen the high-energy main jet exhaust gases- are 
reversed back into the base region. 

This criteria v/.as developed by .studying two-dimensional flow 
across a rearward facing step. Chapman (Ref. 2-1), studying laminar 
flow, was the first to develop a theoretical miodel of the flow reversal' ■ 
mechanismi. Korst, et- al (Ref. 2-2) developed a similar model for 
turbulent flows. '■ Goethert (Ref. 2-3), studying flow reversal from, 
clustered rocket engines, was the first to apply these -theoretical.. models 
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"J" - BOUNDARY OF UNDISTURBED 
,,0ET 

"D" - DISCRIMINATING- LINE, i.e., 
LOWEST ENERGY STREAMLINE 
■ WHICH. CAN- PROCEED- 
DOWNSTREAM 



"D" OUTSIDE "J" NO FLOW REVERSED 



"D" INSIDE "J" . FLOW REVERSED 


FIGURE 2-2. FLOW REVERSAL CONDITIONS' 



cirecciy to the rocket base flow problem. The use of back step flow 
models as diagnostic techniques for predicting the base environment 
has not received wide acceptance because of disagreements obtained 
with CGuivalent experimental results. One of the reasons cited is that 
in order to apply the relatively simple two-dimensional back step flow 
medei to the very complex three-dimensional rocket base flow reversal' 
problem, many significant simpliiications and approximatiop.s mmst be 
made. These are thought to introduce errors in the computed results. 

One of the tasks of this study was to determine the accuracy 
and utility of the flow reversal analytical ’models. For this purpose, 
a computer program was developed and is' discussed in detail in 
Volume III of this report- Comparison between the analytical flow 
reversal model and the- measured results obtained during scale model 
and flight test are shown in Section 6 and in Volume III of this report. 

2. 5 PhUME RADIATION 

The thermal radiation emitted by rocket exhausts depends- on 
.many factors. The basic factors are associated, with the propellant 
corribinations and engine operating conditions, ■ To achieve maximum 
specific impulse (thrust-to-mass-fiow ratio), ‘engine operating co.ndi- 
tions are required which directly influence the plume radiation. For 
liquid-propellant engines, the maximum specific impulse is obtained 
when the mixture ratio is fuel rich. For solid-propellant motors, 
metallic particles are mixed with the fuel to increase the specific impulse. 
Large liquid engines also have very fuel- rich turbopump gases dis- 
charged into the outer periphery of the main exhaust streams near the ■ 
exits. As a result, liquid engines dis.charge large amounts- of excess 
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fuel which can and do mix with the ambient air and* burn under the. 
correct conditions. The afterburning appears to be quite significant 
at low altitudes and produces a high temperature mantle which sur- 
rounds the exhaust plumes. At high altitudes, the afterburning either 
decreases to a very low level or ceases altogether. For liquid engines 
where the fuel is a hydrocarbon- (RP- 1 ), the chemically inefficient com- 
bustion process (i. d. , excess fuel) taking place within the engine system 
produces significant ainounts of particulate carbon (soot). This condi- 
tion, in conjunction with the afterburning, causes the plumes to radiate 
similarly to a high-temperature blackbody. Because of the solid par- 
ticles, the radiation from large, hydrocarbon-fueled engines and solid- 
prcpeilant engines is similar. If the fuel is a nonhydrocarbon Oiydrogcn), 
the exhaust gases are semiopaque and the radiation is, therefore, much' 
lower. 


Rocket exhaust plumes are nonisothermal, nonisobaric, and 
r-onhomogeneous gas bodies. The equation which governs such a 
radiating gaseous body is (see .Figure 2-3) 


qir 


= f.lfs «>«vlbv^-p (-/sP“vdS') ' (2-15) 


where 


q-ir 

P 

I'V 

Ibv 

S, 


• COS 6 dS dco dv 

incident radiant flux 
local density 

spectral absorption coefficient 
Planck function 
optical path 




d - solid angle 

V - wave .number 

6 - see Figure 2-3* 

Cornnuter programs for calculating the rocket exhaust plume radiation 
using Fquation 2-15 are available (see Reference 2-4 for example). 
However, the commuter programs require, . as input, the complete 
spatial ceiinirlon of all. thermLodynamhc and radiation properties along 
all optical paths. Tnis type of information is seldom, available and is 
difficult to generate theoretically for miultiple plumes with afterburning. 

For an isothermal,, isobaric, and homogeneous gas, Equation 
-15 can be reduced to ' • 


where 




nr 


^b ^p 



^dAs - dAp 
P 


(2-16) 


crT^ - plume blackbody ernis.sive power 

f(T"o, p, 1) - effective volumetric gas emissivity 


. dA, 


Cos 0 dco 

7T 


- view factor between dAs and dA^ 


i-Q - plume temperature 


I - mean beam length. 

To .derive Equation 2- 1 6 , the emissivity was assumed to be a function 
of the gas temiperature, density, and some characteristic dimension (i ) 
of the gaseous body. This concept developed by H, C. Hottel is 
discussed in imost engineering heat transfer texts. The characteristic 
dimei'isiozi for circular cylinders and spheres can be obtained from 
I'leference’ 2-5.. The restriction that the gas is isothermal, isobaric, 
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3.nd homogeneous severely •lin'iits thv use o£ this approach when com- 
puting the radiation from rocket’ exhausts • 

Because of the solid particles in' the exhaust plumes, the optical ' 
depth is very short for both hydrocarbon-fueled liquid and solid -propel- 
lant engines operating at low altitudes. Therefore, the radiation- is very 
nearly a surface phenomenon. For this condition, Equation 2- 1 5 'can be 
reduced to the usual equation for radiation between a finite radiating 
surface and a diff erezitial receiver. In finite difference form, 'this 
ecuation can be written as 


^ir 


= E E. 








(2 -.17) 


To use this equation to* compute the radiation, it is necessary to know 
the spatial variation of the plume emissive power (Ep). 


Because of the very complex nature of the problem, accurate 
radiation predictions are very difficult to make. 

2.6 MISCELEANEOUS EFFECTS 

Considering the discussions of the convective 'and radiant heating 
in a rocket base region, some of. the obvious factors which establish 
the base environment are 

o Engine spacing 
o Ease structural arrangement 
o Propellant combination 
o Vehicle operation conditions. 

Tlowever, many other (not-so-obvious) factors also exist- which can 
iiifluence the thermal environment in direct ways. These factors are 
often unique to one specific vehicle configuration and/or to one specific 
flight. Some of these are discussed below. 
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Base heat shields 
iigr-tvveight insulation rr.a 
As a result of the humid 


are usually coated with a high-performance, 
terial. These materials are often quite porous, 
atmosphere in conjuiiction with' surfaces at or 


near cryogenic tem.peratures, a considerable amouart of condensation 
is in contact-with the vehicle prior to launch. It is often possible for 
soTTie of. this condensate to collect above the heat shield and/or be 
absorbed into the porous insulation material. During flight the pressure 
surrounding the base region is constantly decreasing, and thermal 
energy is constantly being added to the materials in the' base region. 

As a result, any moisture changes phase from a liquid to a gas or sub- 
limes directly between the solid and vapor state. The heat shield absorbed 
moisture outgasses directly into the convective boundary layer and 
appears to influence the convective environment during the early phase 
of booster flights. ' Because of the large energy required in changing 
phase, the temperature of the heat shield is lower than that which would 
be expected for a dry heat shield. 


Another example of. phase- change- related environmental effects 
is associated w-ith cryogenically cooled engine walls. Typically, liquid 
propellant' rocket engines are constructed of thin wall tubes brazed 
together along the axial direction. Wall cooling is provided by pum.ping 
some of the fuel through the tubes. Since the Saturn upper stag.es used 
LOX/Hg propellants, hydrogen was pumped through the tubes. It is 
important to the coolant, circuit flow rates that the coolant enter and 
leave the cooling circuit as a liquid. For hydrogen, this means tem.pera- 
tures of approximately -250 to -400® F. Flo.wing hydrogen .at these-tem- 
peratures produces gas side wall temperatures' well below the satura- 
tion temperatures of the mostly water vapor- exhaust products. As a 
result, condensation and even solidification of the water yapor adjacent 
to the wail appears to have taken place. Since- the engine boundary layer 
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fluid is the first portion of the fluid which is;reversed into the base 
region, it would appear that the condensate could have significantly 
influc'nced the temperature of the base gas. 

Extreme engine gimbal patterns and/or inoperative engines can 
also significantly influence the base environment. During a nominal 
flight, the engines girrfbal (i. e. , pivot about an apex) only slightly to 
make minor adjustments in the vehicle trajectory. However, such 
liiings as inoperative engines, gimbal actuatox* failures’ or events which 
create severe vehicle side loads (e.g., loss of the micrometeoroid 
S'hield during launch of the' Skylab) . cause s ever e- engine movements. 

All of these conditions have been experienced during some of the Saturn 
vehicle flights. V/lien these conditions exist, the base environment is 
usually influenced, significantly. As a result, the vehicle base thermal 
protection system design environiment is often dictated by one of these 
more severe 'conditions. 

The existence of soot from hydro carbon propellants is another 
factor which complicates the base heating problem. Since the base 
region is a high radiant heating environment, it would be very desirable 
to use surfaces whiclP have a low absorptivity. Although it is customary 
to employ low absorptivity surface coating, for added insurance, the effect 
can seldom be relied on because of the possibility of degradation result- 
ing from soot deposition. Photographic observations have- shown that 
the base region is almost engulfed in thick exhaust gases, just after 
engine ignition, and occasionally during the release and launch phases 
of the mission. Other photographic observations of the vehicles at 
extreme altitudes have shown what appears to be soot deposition along . 
the cylindrical aitez'body.' This is explained by the reversed gas flowing 
’out of the base region and forward along the vehicle in a highly separated 
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flow regiori. This amount of soot deposition along the side of the vehicle 
irrn^lies an even larger amount within the base region. Needless to say, 
the problern of soot deposition creates serious instrumentation problems 
when attempts are made to measure radiative and convective heating 
parameters. 

Several of these effects will be discussed in more detail in the 
following sections. 
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3. FLIGHT VEHICLE I iVSTRUMEiMTATlON PROBLEMS 

3 . ! INTRODUCTION ' 

Because of. the extreme complexity of the base heating problem, 
one would expect to gain a considerable understanding and insight into 
the problem by studying and comparing the in-flight test results. To a 
degree, the latter is possible; however, certain instrumentation- 
. related problems have thus far prevented the full use of the flight-test 
results. ' , , 

With respect to efforts made in measuring the base environment, 
the Saturn vehicles were probably the 'most well- instrumented vehicles 
ever ilight-tested. Over 200 individual ins-trument types and locations 
were used in the base region. Many of these were carried on several 
of the 30 individual, well- instrumented, flight-test vehicles. However, 
these efforts appear to 'have fallen short of producing -results which can 
be used'to evaluate many of the parameters that are important and neces- 
sary in under standing the base heating problem. If a single factor could 
be identified as a major .cause for the general lack of understanding of the 
base heat transfer, it would most surely have to be the inability to accu- 
rately measure the in-flight base environmental conditions. The latter 
might imply criticism of those responsible .for the design and selection 
of the individual' instrurnents. This is not necessarily intended. The 
truth of the imatter lies in the fact that conditions existing within the 
base region produce requirements which are nearly impossible to meet 
with instruments designed along conventional lines. 

The instrumentation used in the base region was intended to serve 
a dual purpose. First, it was to establish the adequacy of the base ther- 
mal protection system. Second, it was to provide basic information with 
regard to the base heating phenomenon. The types of instruments used 
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wex-c thermocouples, gas temperature probes, pressure probes, a-nd 


n : 
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ment; the high-frequency, high-g conditions; and transient, possibly 
contanninating (i. e. , 'soot from engine exhausts) environments produced 
difficult conditions for the instruments to meet. 

The basic recuir emients for the instrumentatioxi were tliat it 


be lightweight, compact, and rugged. 


:ly in the program it wai 


realize 


;d that accurate, responsive, and reliable instruments were not 


available to meet these requirements. As a result, an extensive instru- 
ment development program was carried on over the entire program at 
significant expense. The instruments used on the later flights represent 
the* most well-developed instruments. However, as will be shown, -even 
these leave much to he desired with respect to individual instrument 
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> . 2 GAS T EMPE RAT U RE P ROB ES 


3. 2. 1 l-'hysical Descriptions 
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Sliown in Figures' 3-1 through 3-7 are sketches of some of the 
pi-obes used to tneasure gas temperatures in the base regions. Gen- 
erally, the probes consisted of, specially designed and constructed 
thermocouples. 
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The probes s.hown in Figures 3-1 and 3-2 were used in the 
Saturn I through Saturn IB base regions. These probes consisted .of 
a chr crnel- alumel thernaocouple mounted in an insulating case. 

Various radiation sbields were used in an effort to reduce the effects 
of absorbed radiation on the temperature readino;. On the Saturn IB 
vehicle, the probes were mounted so that the' distance from the surface 
to the probe sensing element was varied as shown in Figure 3-3. 
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'IGURE 3-4, SCHEMATIC Or SATURN S-IC ENGINE GAS TEMPERATURE PROBE 
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Shown in Figures 3-4 through 3-6 are sketches of the' probes 
used in thie Saturn V, S-IC base region. The probe shown in Figure 
3-4 was used exclusively aro‘und the engines, while that shown in 
Figure 3-5 \vas used on the heat shield. The heat shield probes were 
•mounted trirough the heat shield with the distance from the surface to 
the sensing element varying between 0.25, 1.0, and 2.5 inches. 

As seen in the figures the designs of the sensing elements of 
the S-IC heat shield and engine gas temperature probes are distinctly 
different. The engine probe consisted of a bare 0. 002-inch-diameter 
platirrum/lO percent rhodium thermiOcouple surrounded by a double 
radiation shield. The inner and outer radiation shields were plasma 
spray coated, with aluminum oxide to prevent catalytic reactions. Gas 
. contacted the thermocouple junction by traveling a circuitous rovlte 
throvigh a 0. 15-inch square hole in the outer shield and through 0.040- 
inch or 0. 10-inch-diameter holes in the inner shield. 

The heat shield probe consisted of a fine platinum/ 10 percent 
rhodium thernaocouple v/ire encased in a platinum/6 percent rhodium 
sheath filled with compacted magnesium oxide insulation .(see Figures 
-5 and 3-6). To reduce thermal radiation, the sensing element was- 
surrounded by a single outer shield. To prevent catalytic reactions, 
the exposed surfaces were coated with molybdenum silicide, and the 
sensing element was coated with ceramic cement. Gas entered the 
cuter shield through three equally spaced radial holes 0. 1625 inch in 
dia-meter (approximate). 

Figure 3-7 shows a schematic of the gas temperature probe 
used in the base region of the Saturn S-II stage. These probes were 
unshielded and extended approximately 2 inches from the heat shield 
surface. The design of the sensing element which contained the 
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thern'ioccupie was considered proprietary by the manufacturer, 
ever, the information which could be obtained suggests that the design 
of the S-Ii staere sensing element was similar to that of the S-IC heat 
shield probe showm in Figure 3-6. No doubt different materials such a.s 
the aluimina sheath were used in construction; otherwise, the design was 
probably sinvllar. 

Further information regarding details of constructing and quali- 
fication tests o’£ the Saturn gas temperature probes can be obtained froim 
References 3-1 through 3-8. 

3.2. 2' Gas Temperature Probe A^ccuracy 

Based upon the data available, it must be concluded that the 
ability of the gas temperature probes to measure accurately the tempera- 
ture of the base gases is very questionable. For some conditions the 
probes may have indicated the true temperature within a few percent, 
however, this is a result of fortuitous circumstances. Furthermore, 
it does not appear possible at this time to correct the erroneous inflight- 
measured results with any degree of confidence and/or accuracy. The 
reasons for these unfortunate conclusions are many and varied and are 
discussed in detail in the following sections. 

Briefly stated, the gas temperature probe errors appear to be a 
result of the following factors.: 

’o Coating the surfaces of the probes to prevent catalytic 
reactions produced surfaces with high radiation absorb- 
tivities and caused the probes to be overly sensitive to 
Incident radiation and reradiation at high temperatures. 

o Capsulatihg the thermocouple inside the insulated sheath - • 
would appear to have had sevez'al deleterious effects: 
the heat transfer path be'oveen the gas and the thermo- 
couple was increased; the sensible heat capacity of the 
therrnocouple/insulation/sheath system was considerably 
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higher than it would have been for the thermocouple 
alone; and the conduction path beuveen the body of the 
probe and the theriTiO couple was increased. 

The net effect of these factors could have caused the propes to measure 
both higher or lower depending entirely upon the instantaneous equilib- 
rium conditions which prevailed. V/Tien these probes are used in the 
base region some, of the effects' tend to cancel each other, providing 
some reason to speculate that the overall accuracy might be better than 
that indicated by considering each source of error separately. 


^ Gas Temperature Probe Theoretical Eval 


uation 


Considering the schematic of the gas temperature probe shown 
in Figure 3-8, \ a heat balance about the thermocouple junction gives 
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Qg - net heat in the junctioi 


- ' net radiation absorbed or: rejected by the junction 
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- convective heating to the- junction 


heat conducted into or out of the junction. 


j-\. s s u m iTi.^ u 


the junction to be an isothermal mass, the net heat in the 


junction at any ti'me is the sensible heat; therefore- 


Qs = ""j cpj -3; 


(3-2) 
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v/j - weight of the junction 


C-o: - specific heat of the junction 
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IGURE 3-8. SCKhMATIC OF A TYPICAL THERMOCOUPLE GAS TEMPERATURE PROBE 
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Tj - tornpcraturc of the junction 
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Considex'ing the radiation interchange betv.'een soine small but 
finire junction area, xlA.-, and an arbitrarily picked surrounding area, 


A A,, , th'’ 


xe radiant energy emitted by AAj which is absorbed by AA^. is 


Qj-K 
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A A; Z 


j ^b; 


(3-3) 


v/ner e 


Q; 


- net radiant energy emitted by 'AAj and absorbed 'by 


A A. 


- radiation interchange factor 




- biackbody emnssive power of AA: 


J 


ft-, 


= cr T 


•oj - -J • 

The radiation interchange factor includes not .only the direct radiation 
belnveen AAj and AAj^ but, also, any reflected and rereflected energy. 
The net radiaxit energy exchange,’ i. e^ ,■ the amount emitted by AA^ and. 
received by AAj niihus the amiount received by A Aj^ that was emhtted 
by AAj (Ref. 3-9), is given by 




(3-4) 


Then, if ail of the surrounding (which at this point includes the shield, 
radiant sources, etc. ) are divided into n arbitrary area elements, the 
net exxergy exchange between AAj and the n radiating elements is 
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(3-7) 


Then, using the reciprocity theorem (Ref. 3-9), 

3^K-j “ '^j-K (3-S) 

gives . , 

?'j-K - %j) • ' (2-9) 

Letting the thermocouple junction he subdivided into m elemental 
areas, AAj, the net radiant energy exchange between the junction 
and all of the surroundings is 

m n 

Qr = Z Z Pj-K- - =bj> • (3-10) 

j = lk' = i 

Equation- 3-10 can be simplified by the following: 
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Then, s.ince 
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the second terra in Eauation 3-Ii can be written as 
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Then, assuming £j E-j^j = £j c*Tj to be constaiit over the junction, Equa- 
tion 3-13 beco-mes 


m 


) AA; £ ■ Ev - = ■ £ • Ev • y' A A . 

.) ^3 hj ■ .^3 -^&3 L .‘-^3 
3=1 3=1 


e 


= A; £; 0- T; 


•j ^Oj - C3 


(3-14) 


{ .'! 


n 

L * ' 


U 

D 


If 1 

u 

( 1 


whe 


re 


A . 




- surtace area oi the junction 
-■ emissivity of the j'onction 


Tj - temperature of the junction. 

Equation 3-14 is the thermocouple radiant, emission ternr and 
accomits for all of the energy lost by radiation that is absorbed by the 
surroundings. Thus^ the net radiation can be written as 
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T'iC first ternr on the ri<>>it of Ecuation 3-15 re*Dresents all the 
which radiate to the junction. For convenience, this can be subdivided 
into the separate sources of radiation. Therefore, dividing the plumes 
into p radiant sources, .the radiation shield into f sources, and letting a 
represent all other sources which radiate to the junction gives 
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The convective heat flux to the junction can be written as 


O = A- h* - T-) 


(3-17) 


wher e 


= r T 


r - recovery factor 


Then, substitution ofEquations 3-2, 3-16, and 3-17 into Equation 3-1 
gives 
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sensible heat capacity of the junction 


•CiT: 

'"j dr 


= V ( y aa.at. b 

mL-' • 

. j=i\g=,i 


■ S. ^j-S 

0=1 e=l 


e -^he 




The indivi-duai terms of Equation 3-iS reoresent the following: 
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net radiation ei'nitted by the plumes 
that is absorbed by the junction 
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net radiation emitted Dy the radia- 
tion shield that is absorbed by the 
junction 


'( i I 


E E/ 

1=1 e^i 


i -r-' Tn- 

*j'^''j-e ^be 


A . ^ . 0“ 1 • 
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net radiation emitted by all other 
radiant sources that is absorbed 
by the junction 

net radiation given up by the junction 


convective heating to the jvuiction 


conduction heat lost or gained by 
the junction . 


n- ! 

U I 


c \ : 

: 
i 


d_j 1 


U I 


For an ideal probe, all terms of Equation 3-18 are identically 
zero, giving T; = i* T^y « T ; thus, measuring the emf output of the 

Jag 

qrobe gives the approximate gas temperature. For an actual probe. 
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all terrr.s in Zouation 3-18 are not zero and the result is 
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likewise, 


T. ^ r T, 


By studying Equation 3-18 it can he seen that in the design of an 
accurate gas temperature probe it is desirable to have the following: 

o The sensible heat capacity of the junction (or -specifi- 
cally the mass of the junction) should be as small as 
possible [ i. e. , wj Cpj (dTj)/(dr) 0] . This gives a 
rapid response to changes in temperature. ’ 

o The conduction path between the therimocouple junc- 
tion and the iprobe. body should be as small as possible. 

This reduces the heat lost. or gained by the junction via 
■ conduction (i. e. , Q,^ ^^0). 

G The-net effects of absorbed and emitted radiation should 
be made negligible. 

o The gas should contact the probe junction v/ith low flow 
velocities so that the recovery factor, r, would be as 
close to unity as possible. , 

o The heat transfer coefficient should be made as high 
as possible so that the convective flu>: to the junction 
(and therefore the sensitivity to the gas temperature), 
would be high relative to the other sources' of. heat 
transfer. 

Trie requirements for which, the Saturn gas probes were designed, such 
as the ability to withstand high-g, high-frequency vibrations; high- 
absorbtivity surface coatings-; etc., made these idealistic requirements 
almost inroossible to an’oroach. 
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Shielcl'hig the junction reduces the net effects of absorbed and 
emitted radiation. 3y following the derivation given on page 311 of 
Reference 3-'9 for the'-net radiant flux transferred between two shielded 
surfaces, it can be shov/n that the net flux transferred between the outer 
shield and the inner shield is given by 




^net 


.y ^ 

k “^1 


(3-20) 


‘on 


where 

n - number of surfaces between the junction and the outer' 
shield including both. of these .surfaces 

emissive power of the outer shield 

p - reflectivity 

£ - emissivity. 

Considering a single- and a double- shielded thermocouple (see 
Figure 3-9), the effect of shielding the junction can be shown. For the- 
single- shielded thermocouple junction, Equation 3-20 gives 
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For the double- shielded jionction, Equation 3-20 gives 
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It can be seen from Equations 3-21 and 3-22 that even if the sur- 
face reflectivities were zero (i. e. , biackbody surface), using a double 
shield reduces the net flux by a factor of two below that of the single- 
shielded probe- Experimentally, the effects of radiation- shielded 
thermocouples and- conduction influences are reported on pages 263 and 
264 of Reference 3-10- 

3.2.4 Gas Temperature Probe Qualification Test Results 

Qualification tests conducted on the Saturn gas temperature 
probes (Refs. 3-3, 3-4, and 3-7) show conclusively that the in-flight 
rnea.sured temperatures must be considered questionable- The purpose 
of the aualification tests was to demonstrate by actual testing that the 
instruments could meet the mechanical and functional requirements- 
The thermal tests consisted of exposing the. probes to a hot gas flow 
and a radiant flux. The orientation of the probes as v/ell as the flow 
velocity, gas temperature, density, and radiant flux were systematic.ally 
varied. 

A typical hot-gas test consisted of exposing a heat shield probe 
to a flowing gas with known velocity, stagnation temperature, and 
density. A test fixture was designed so that the probe could be rapidly 
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injecrcd into- the gas stream once the gas flow condition was established. 
The temperature of the probe was monitored until well after steady- 




Some tynical results obtained during the qualification tests arc 
shown in Figure 3-10. After exposure, the probe would z^ise to a 
steady- state value that was well below the actual gas temperature.. 

For condition 4 of Figure 3-10, tiie actual gas temperature was approxi- 
mately 2054'' F; at steady state the probe indicated a temperature, of 
llOO'^F. During this test the flow velocity was .1928 ft/sec, the stag- 
nation pres suzm was O'. 81 psia, and the density ’was 0. 00097 Ib/ft^.' 

The S64‘''F temperature difference represents a probe error of approxi- 
mately 34 percent (based on degrees Rankine).. ’ ■ As may be seen in the 
figure, the errors are less at higher flow rates and lower temperatures. 

In addition to varying the mass flow rates, the probes were 
injected iiito- the flow with different orientations with respect to the 
centerline flow direction. Probe positions were. varied from: probe 
axis normal to the flow (Position I); probe axis 60 degrees to the flow 
(Position II); probe'axis 30 degrees to the flow (Position III); and probe 
axis parallel to the flow (Position IV). 

The results of approximately 44 tests conducted in the hot-gas 
• cnwlroizment are reported ^in Reference 3-4. These test results were 
evaluated to detezmnine the probe error as a function of mass flow rate, 
position, azid gas remmerature. The test results are tabulated in Table 
3-R The probe temperature, Tp, is the apparent steady- state value. 
A-lso shown in Table 3-1 are the average conditions indicated for all 
Position I tests. These results are plotted in Figure 3-il. Across 
the abscissa of this- plot are- shown the approximate times that the 
flow in the S-IC base region had the same mass flow flux (pu). 
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FIGURE 3-10. TYPICAL HOT GAS FLOW. TEMPERATURE PROBE QUALIFICATION TEST RESULTS 
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The test results shown in Figure -3- 1 1 were obtained by exposing 
tiee probe to hot gas only. The test chamber was a water-cooled pipe. 

The temperature of the chamber wall is not reported in Reference 3-4; 
however, some radiation exchange between the pipe wall and the probe 
could have influenced the results. 

During another part of the qualification tests the probe was 
exposed to a radiant flux in the absence 'of any forced convection. The 
results of these tests are shown in Figure 3-12. When the probe was 
irradiated with a 40~3tu/ft“ -sec flux from all sides the probe indicated, 
a steady-state temperature of 2500°F. Exposure to a 15-Btu/ft’ -sec 
flux (approximately the magnitude of the S-IC heat shield incident flux) 
produced an. indicated steady- state temperature of approximately MOO^'F. 
W'hcn the sides of the probe were shielded with a special reflective 
shield and the end was exposed to a 40-Btu/ft’ -sec radiant flux, 'the probe 
indicated a steady-state temperature of approximately .450'' F. 


A similar set of tests was conducted using the S-IC engine gas 
temperature probe. The engine probe was exposed to hot gas at vary- 
ing ffov/ rates and to combined radiant fTox and gas flow. During the 
hot-gas flow tests the gas temperature was varied between. 300 and 2100"* F 
(Figure 3-13). During the combined radiant and gas flow tests the probes 
were exposed to incident Ineating rates of 19.3 and 10.6 Btu/ft^-sec in a 
gas flow having temperatures of 20 to 60'' F. The probe was also exposed 
to a 400®F gas tetmper ature flux in conjunction with the 1 0. 6-Btu/ft^ -sec 
radiant flux. 


Based upon the' steady-state conditions, the engine probes con- 
sistently indicated low when exposed to hot gas flow only, and high when ' 
exposed to the combined radiant and gas flow. During the radiation tests, 
when thc 2 :>robe was exposed to a 10. 6-Btu/ft^ -sec radiant flux, the indi- 
cated error was less when the gas temperature was 400" F than when the 
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gas ten'pcrature was between 20 and 60T. The latter rnay appear 
to be somewhat misleading since it indicates improved accuracy, at 
iiigher gas temperatures; however, these results actually indicate that 
with the probe exposed to a 10. 6 Btu/£t“-sec radiant flux the junction 
equilibrium teimoerature is closer to the 400° F gas temoerature than 
to the 20 or 60° F gas temperature- This also illustrates that the probe 
errors tend to be reduced when the probe' is exposed to a combined radi- 
ant and het-gas environment. 


The Saturn S-II stage gas temperature probes (Figure 3-14) were 
also tested in a gaseous environment (Ref. 3-7). These tests coa-isisted 
of exposing the S-II gas temperature probes to the flame of a plasma, 
spray nozzle. The actual gas temperature was measured with a so- 
called standard thermocouple located in close proximity to the probe 
juncricn. All tests were. .conducted at sea-level, ambient pressure, 
and the results are shov/n in Table 3-2. As can be seen in the table, 
the temperature- indicated by the probe is significantly les s than that • 
indicated by the standard, particularly at the high-gas temperatures. 
Another factor which should be considered in these tests is the prob- 
able increase in error which w.ould occur with decreased density. Hot 
gas tests of the other' Saturn probes indicated an increase in error as- 
the mass flux (p'u) was decreased. Since the S-II stage probe tests 
were carried out at sea-level, ambient pressure, it .s probable ..that • 
the density and likev/ise, the mass flux (pu), were considerably higher' 
than that which would have been associated with the S-II stage during 
flight* Therefore,, the error in the gas temperature measured by these 
probes during the flight might be higher than that suggested by Table 
3-2. . • • 
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TABLE 


3-2. SATURN S-II GAS TEMPERATURE PROBE 
QUALIFICATION TEST RESULTS 


GAS TEMPERATURE INDICATED 
BY STANDARD (=F) ■ 

GAS TEMPERATURE INDICATED 
BY PROBE (°F) 

DEVIATION 
( = F) 

1072 

930 

• 142 

.1333 

1210 ■. 

133 

1821 

1530 

291 

217.2 

1 800 

■' 372 

2879 

2340 

539 


During this study attempts were’ made to correct the flight data 
by developing very detailed thermal models of the gas • temiperatur e 
probes and surrounding heat shield. The measured flight test data was 
intended to be used as.iuputs to the analytical model, ' These attempts 
were abandoned' when, in the course of the- study, it was concluded that 
so many of the variables necessary for a solution could not be specified 
with any degree of accuracy and that the computed gas temperature would 
be at least as questionable as the measured value. As an example, the 
heat transfer coefficient for the heat shield probe sensing element could 
be evaluated fro'm the qualification test data. 

Figure 3-15 shov/s a dimensionless Nusselt number/Reynolds 
number plot of the convective heat transfer to the probe sensing element 
obtained from the qualification test results. Also shown in this figure 
is the convective heating to a cylinder normal to the flow obtained from 
Reference 3-10. As seen in the figure, the qualification test results indi- 
cate (approximately) an order of magnitude variation in the convective 
heat transfer coefficient. The convective heat transfer to a cylinder 
normal to the flow is one to two orders of magnitude above the results- ^ 
obtained from the qualification tests. The large variation in convective 
heat transfer to the probe sensing element was possibly caus'ed by the 
radiation shield. . The large deviation in the convective heat transfer 
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coefficient- v/hich would result from this data would make the computed 
results very questionable. 

Sudunas (Refs. 3-5 and 3-6) evaluated the S-il. gas temperature 
probe, results using an analytical model of the probe. In this analysis 
t.he probe convective coefficients were computed by a circuitous route 
v/hich involved relating the convective fl'ox measured by the heat flux 
meters to the probe coxrvective flux. Several important assumptions 
were necessary in making the analysis and the final results depend upon 
the ccrrcctness of these assunxptions; however, the computed gas tem- 
peratures were within :i:100°F of the measured gas temperatur es . The 
analysis of the S-II stage probes was considerably easier than that 
z'equired for the S-IC probes because of the absence of a radiation shield 
and the low incidentmadiation in the S-II base region. 


3.3 - - TOTAL HEAT FLUX METERS 
• 3-3.1 Physical Description 

During the course of the Saturn flight test programs several 
different types of total and radiant heat flux meters were used. Sketches 
of some of the differ ezit meters and techniques used in evaluating the 
heat flux from these meters are contained in Volume II, Section 4 of this 
report. Briefly, two basic types of heat flux meters were used to mea- 
sure the total heat flux in the Saturn base regions. Some of the charac-' 
tcristics of these meters are presented in the following paragraphs. 

Early in the program the so-called slug-type heat flux meter 
was used. This instrument consisted of a thin, rnetai slug which was 
insulated from its* surroundings and a thermocouple ^attached to the slug. 
The heat flux was evaluated by using various calibration methods which 
were related to the slope of the measured temperature-time history. 
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:/f tr.cter was used early iu the Saturn program* Continued use 


)scc scv'crai or o ole ms tnas v. 


inherent to the slu^x-tvpe metcr. 


T?:c slug could not be completely isolated from the mounting structure; 
thus, the meters had sizeable conduction losses, particularly at high 
temperatures. Changes in the thermal properties of the surface coat- 
ings limited the rnaxiimum teimperatur e to which the meters could be 
calibrated. Often, during flight use, the mieter temperature would 
exceed the m;aximium calibration temperature by several hundred 
degrees, requiring extrapolation of the calibration curves. This caused 
somic doubt as to the validity of the extrapolations because of possible 
changes in the surface properties and the higher conduction losses. 

Also, because of the tcmperatui^’e limitations' and the high-g environ- 
nnent, a relatively large slug was used. This caused the instrument 
to have a relatively slow response to the highly transient base 
c n V i i' o n n n cut. 
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Later in the Saturn program an instrument was developed which 
did not empear to have most of the undesirable ciualities of the slug-type 
meters. This instrument used, a completely different principle for 


measuring and evaluating the hedt flux* This instrument, schematically 
shown in Figure 3-16, consisted of a thin, constantan foil disk bonded 
at its periphery to a copper cavity. The copper- constantan formed a 
thermocouple pair which generated an emf in proportion to the tempera- 
ture, difference between the center of the-foil and the edge. The tempera- 
ture difference is almost proportional to the heat flux incident upon the • 
,foil. The naass and heat capacity of the copp’er was made relatively large 


in compai 


risen to the mass and exposed 


of the foil. As a result. 


the heat absorbed by the foil was rapidly conducted into^.the copper. Even 
w’lien- exposed to a relatively high heat flux, this technique prevents the 
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TYPICAL SATURN S-IC TOTAL HEAT FLUX METER CHARACTERISTICS 
RANGE; 0 to 60 Btu/ft^-sec 

0.06A5 in. < B < 0.0655 in.. 

S n/ 0.003 in. 

r- • 

Q = Z) Ap c'^ -> (Btu/ft^-sec), C - emf output 
n =0 

Ao = -0.102, Ai .= 6.462, Aa'' = -0.00959 

A 3 = 0.0-0149, ki, =-0.0001 16, As = 0.358 x 10 


FIGURE 3-16.- SCHEMATIC OF THIN-FOIL HEAT FLUX METER 


cci'itcr terriperatui'e of the foil from rising more than one or two 
hundred degrees above the temperature of. the copper sink, which 
changes only slightly. ' - 

This t>me of instrument can be made, very responsive to chang.es' 
in hear fJ.uM. (The Saturn instruments had. response times on the order 
of 0, 150 second. ) Another advantage-of this instrument is that the heat 
flux is almost linear with emf over a wide range of heat flux. These 
instruments were also calibrated and a curve fit of the results has 
shewn that the incident flux for instruments of the same type could be 
•obtained using a fifth order polynomial with variations of less than 
±0. 3 percent. 


Because the surface temperature is important in evaluating the 
convective flux heat transfer coefficient and cannot be determined from 
the caloriimeter emf output, a separate thermocouple was attached 
between the foil and the copper (Figure 3-16). 


3.3.2 Total Heat Flux Meter Accuracies , 

¥/iien slug-type heat flux imeters were utilized, three different 
techniques were often used to evaluate the heat flux. Each of these 
techniques gave different results which should be compared when con- 
sidering the ‘accuracy of the imeters. Several comparisons of the 
measured flux evaluated using the different methods are shown in 
Volume II (see Figure 7-21 for an example). The data band' represented 
by these methods has often been taken as the realistic accuracy of the 
individual instrunrents . However, based upon the techmiques used, it- 
would appear that the preflight calibration results should be more 
realistic during the early portions of the flight, and the cutoff decay 
results more reliable- during the latter phases of the flight. 
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The stated accuracy of the thin-ioil total heat flux meters is 
±10 percer.t; however, the qualification test results (Refs. 4-16 through 
4-10 of Volume II) sho\v that accuracies considerably better than ±10 
percent can n.crrnally be expected fromi these instruments. 

3.3.3 Total Pleat Flux Meter Surface Temperature 
- Mismatch Problem 

One factor which may have affected all of the total heating rate 
measurements is the so-called ^heimperatur e mismatch problem’^ This 
■orobiem, shown s chematically in Figure 3-17, is a result of the tempera- 
ture differ ence-.between .the surrounding insulation and'Mhe heat flux sur- 
face.- This problem, is that if,- during the flight, the' surface teimperature 
of trie insulation is significantly different from the surface temperature 
of the heat flux meter, the gas boundary layer is perturbated as it 
crosses the meter; and, since the convective flux is related to the 
temperature gradient .and the thickness of the' boundax-y layer,’ it is also 
norturbated. The available experimental and theoreticaT investigation 
into this problem (Ref. 3-11) indicates that the Saturn total, calorimeters 
probably indicated a higher heat flux than would have actually existed 
if the surface had been isothermal and at the same temperature as the 
heat flux meters. Experimental support .fbr. this conclusion can be seen 
in Figure 3-18.’ ‘ 

Figure 3-18 shows some of the experimental results reported in 
Reference 3-11. These tests were conducted by mounting a membrane 
calorimeter in an isothermal surface and using an oxyacetylene flame 
to ccnvectively heat both .the isothermal .surface and the calorimeter. ■ 
The location of the calorimeter from the leading edge wa^ varied from 
6 to 31.6 inches. Shown in the figure is the ratio of calorimeter heat 
transfer coefficient to the undisturbed surface- coefficient. The 
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The reference convective flux (i, e. , convective flux to an LUidisturbed 
curiaco) is shewn along t'ne abscissa. The convective fho:': coefficient 
indicated by t'he naerrfDrane calorinreter was approximately 1*2 to 1.5 
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Some of the Saturn night test data also tends to support the 
•possibility of higher convective flux measvired by the rhennbrane calorim- 
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the Saturn I, Block II flights, one of the instruments 


was changed from a srug-type heat flux meter to a membrane-type 
meter, resulting in a different surface temperature. Figure' 3- 19 shows 
the predicted surface temperature of the insulation and the measured 
tcm-pcraturcs .of the heat flux meters. ■ According to the temperature • 
mismatch theory, the heat fluo: indicated by the membr ane-type caloriim- 
cter snouid iiave been higher than that indicated by the slug-type 
calorimeter. Figure 3-2 0 shows the total heat flux measured by thes'e 
meters and tends to indicate that the membrane-type meter did give a 
iiighor indication. 

3.3.4 Hea.t Shield Moisture Effects 


The effects of moisture absorbed by the heat shield was another 
•factor which probably affected the total heat flux meter measured results, 
Prior to flight moisture could be and probably is absorbed, by the porous, 
heat shield material. During flight the decreasing pressure and heat 
■transferred to the heat shield vaporizes the moisture. The moisture' 
leaving the heat shield insulation e*iUers the boundary layer and contacts 
the total heat flux meter. Since the temperature of the membrane- type 
calorimeter is below the saturation temperature of the moisture (see 
Figure 3-19), the possibility exists that the vapor will condense. Con- 
densation of the moisture on the calorimieter surface would have a pro- 
nounced effect upon the* rrfeasured heat flux. This problem was investi- 
gated for NASA by The Boeing Company and the results are reported 
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During the Saturn progran-: the radiant h 
these instruments in a manner similar to 
iux meters. 

3.4.2 \V in d o w T r an s rr. i s s i o n D z' r o r s 


eat flux was evalviatcd 
that used for tlie total 


During this investigation, a previously undetected error in the 
measured fJ.ight radiant flux was discovered- This error is related to 
the nzanner in v/hich the radiation calorimeters are calibrated, the 
angular dependence cf the sapphire window, and the way in which the 
instruments are used on the vehicle. Theoretically, it can be shown 
that the radiation indicated by the S-IC radiation calorimeters is at 
least iO percent higher than the actual incident radiation. An explana- 
tion of the reason for introduction of the error is presented in the follow- 


ing paragraphs. • ' 

During calibration a standard water-cooled calorimeter and a 
test instrument are simultaneously exposed to an electrically heated 
graphite plate. The graphite heater is specially designed to radiate 
spectrally like a very good greybody- Both instruiments are placed in 
close proximity to the. plate, giving a view factor cf very near unity 
(i. e, , the view factor from the front surface of the test instrument is 
unity, net the view factor to the' sensing element). 
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Because of the window transmissivity and its angular dependence, 
the radiant energy reaching the sensing ele‘ment is much-iess than that 
reaching the standard calorimeter. Calibration of the instrument is 
carried out by equating the emf output of the test instrument to the heat 
fivuc being measured by the standard calorimeter. This equates the 
radiant energy which was absorbed by the sensing element to that which 
was incident upon the front of the instrument. 
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If the radiation which is excluded from the radiation sensing 
element is considered in some detail, it becomes obvious, how the 
errors occur. First, the spectral transmis si\dty of the sapphire 
window (in the normal direction) is high (t ~0.,9) between \ 0. 2 

micrometer and micrometers. Beyond these wavelengths the 

tr ansrnissivity rapidly drops to zero; therefore,- any radiant energy 
emitted by the calibration source which is less than X ss- 0, 2 micrometer, 
greater than X 5. 0 mher ometer s, and approxiiTiately iO percent of that 
between these limnlts is excluded from the sensing element by the window. 
Second, the sensing element is slightly below the front face of the meter 
and, therefore, viev'/s the calibration source through a hole. The result- 
ing effective view angle through which the sensor sees the calibration 
source for the S-IC calorimieter s was approximately 150 degrees. There- 
fore, any radiant energy approaching the calorimeter at an- angle of more 
than 75 degrees to the normal would not reach the sensor but would be 
intercepted by the calorimeter housing. Third, the transmis sivity of the 
'sapphire window norm.al to the surface is appi’oximately t ~0.9; how- 
ever, window transnris sivity has an angular dependence. When the angle 
to the normal is 60 degrees, the transmissivity had dropped from r ~ 0. 9 
to T '^ 0. 775 and at 70 degrees, r "^.0. 6. When the angular transmissivity 
is integrated over the entire effective view angle of a radiation calorim- 
eter the effective transmissivity obtained in r '^0.79. Therefore, the 
difference between r 0. 9 and r 0. 79 is an additional amount of radiant 
energy which does not reach the sensing element. 

When the radiation calorimieter s are used to mieasure the radia- 
tion in an unknown environment, the three excluded radiant energy 
sources are automiatically included in the calibration. If the view 
factor to the radiation source is near unity, and the source radiates 
like a greybody, then the incident radiation indicated by the instrument 
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is c.nproMirr,atciy correct. However, if the instrument is exposed to a 
source which has a view factor significantly less than u:vLty, ai'i error is 
introduced by the calibration in which excluded radiant energy had been 
automatically included. The Saturn heat shield radiation calorimeters 
experienced’ this probioni, resulting in an error which has been esti- 
mated tc be approximately 10 to 12 percent. 


An additional error can be introduced if the radiant source does 
not* emit like a greybody, Greybody radiation suggests that the source 
emits radiation which is spectrally a constant percent of the cori*espond~ 
ing blackbody radiation (i. e., the emissivity is independent of wavelength). 
The calorimeter window- allows only radiation between 0. 2 micrometer 
/GAv A 0 micrometers to pass through; however, the excluded energy 
is included in the calibration because radiation ove;r all wavelengths 
is absorbed by the standard instrument. Therefore,, if the test instim- 
ment is exposed to a radiant source which emits radiation significantly 
different from greybody radiation in the regions 0 . 2 . micrometer > 

X <.5. 0 micromieters (i. e., outside of the spectral range of sapphire), 
tills will not be accurately indicated in the calibration. For vehicles 
using hOX/RP-i propellants, this may or may not be a significant' 
source of error since Aie plumes appear to emit approximately grey- 
body radiation (see Figure 6- 13b). However, for vehicles using LOX/ 

H 9 propellants, the error could be significant because the radiant 
energy from these plumes is not of the greybody type (see Figure 6-13a). 












3 . 4- 3 Derivation of the Window Angular Transimis sion Effects . 

Figure 3-23 shows schematically a calibration test arrangement 
thac depicts the window angular transmis sivity dependence* The radiant 
energy emitted by dA. that passes through the window and is directly 
incident upon dA> is given by 
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.FIGURE. 3-23, SCHE?4^TIC OF RADIATION CALORIMETER CALIBRATION . j 
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-A, = 'i'v/(9 2» >^) (?l > >^ ) . 


(3-23) 
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- radiant energy emitted by dAr and incident upon 


dA-, . 
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window transmissivity 
radiant intensity of dA^ 


^ solid angle intercepted by dA, as seen from dA^ 

Using L-ambert‘s cosine, law 


(6l , \ ) = I, (X) cos <!)1 


the Stephan-BoltzinLann relation 
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and the definition of solid anomie 
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Iquation 3-23 can be wz-itten as 


I I 


'■'w (o 2 > M cos b, cos 6, dA^ dA^ dX 

dQ_A^ _A, ~ __ _ 2 * (3-26) 




3y making the following substitutions 
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arid (sec Figure 3-23) 


G-, =09 = o, = o. 


cA, = Z- de-db , 


COS'* o 


arid 






cos" o 


vouation 3-26 becomes 


*^A, -dA, 


f. (4>, X) E-j_ (X) cos 9 sincp d? d0 dX 

"X ‘"e *^4) , ^ 

(3-27) 
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q - radiant energy emitted by A, that is directly 
incident upon sensor differential element dA 2 


0 


- angle about the axis..- (see Figure 3-23). 


Fouation 3-2? is the basic- equation giving the radiant energy emitted 
by the. calibration source that passes through the window and strikes 
the sensing element’. 

For the standard instrument which is assumed to be a minor 
image of the test instrument with = 1, Equation 3-27 can be inte- 
grated in closed form to give ■ . 
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X=0 


■^1 (X) clX - Cl cr 
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(3-2S) 


, .Equaticn 3-27 for the radiant flux meter cannot he separated 
because the -transmissivity of the window is a function of the polar- 
angle, 6 , and the wavelength, \ . Shown in Figure 3-24 are the sensor 
absorptivity, and the window transmissivity and refractive index.. 
From this figure it can be seen that the transmissivity -can be roughly 
approximated by a square wave with 


0 ^ . X ^ 0 . 2 u , 
0.2p~ X = 5.0 ijl, 
5*. 0 ^ X = CO, 


'^w' = 0 . ■ 


O.S4^r^5 0.9 


(3-29) 
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Since the transniissivity of the window is a function of the polar • 
c^nglo, o , this must also be considered. Figure 3-25 shows how an inci- 
dent beam of radiation is' refracted and reflected as it passes through the 
window. A portion of the energy is reflected at both the front and rear 
s'urfaces. of the window. 

From the theory of optics (Ref. 3-16), the angular transmittance 
is given by 


T (g ) = T ' (6 ) [ i - p (<? ) ] ' 


(3-30) 
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9(6) = ‘T 


sirr ( 6 - 6 ^) 
sin" {6 -r 6^) 


^ (6 -o') 

tan“ (o - 9 * ) J 
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N T i 
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(I - Pr.yj 


(3-33) 


(3-34) 


(3-35) 


t(6) - angular transmittance 

- transmittanc e in the normal direction 

- angular transmittance without reflections 

p (o) - angular reflection • ; 

- normal reflection 

N - refractive index 

K - absorption coefficient 

X - thiclcness in the normal direction. 

Using a window thickness of 0. 020 inch and normal transmiss ivi- 
ties of = 0. 84 and = 0. 9, the radiation calorimeter window angular 
transmissivity v/as computed using Equations 3-30 through 3-35; the 
results are shov/n in Figure 3-26- 

Equation 3-27, which gives the radiant energy incident upon the 
sensor, still cannot be integrated directly. Letting r an effective 

e^i- 

transmissivity of the window which depends on the geometry of the source 
with respect to the test instrument. Equation 3-27 can be written as 
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A- 3 . U pQ-^~ p 92 

E^ i\) d\ \ — - 'I T ( 6 ) COSO sin dtp 

\ = 0 . 2 .^ ^^8 = 0 ' 9 i . '. 


0 ^9— 2tt i^o-> , • . ' • 1 r\ 

( t- /, s -,. f ( cos o sin o cio d0 

'e££ \ \ \ :: 

'"X=0.2 ‘^ = 0 o, 


p\=5.0 

= ^et'£ ^A.-dA, \ Si (X) dX- . (3-36) 

•\=0.2 ■■ ,- 


Then by comparing the first and third equations and integrating, over 0 
it can be seen that 


Oo 


t( 6 ) cos 6 sin cp <ip " 
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‘ eii 


Ai -dA2 


(3-37) 


where is the effective transmissivity of the window considering 

angular dependence.’ ' ■ 

Numerically integrating Equation 3-37 over the effective view 
angle of the calorimeter gives 


■ (d = 90) 

= 0.90, Teff = 0. 7895 

(p = 90) 

T-, = ' 0. 84 , T - = 0. 734 . 

“ eii 

From the results above it can be seen that the transmissivity of the 
calorimeter window is reduced significantly as a result of the angular 


dependence. 

During calibration the standard instrument and the test calori- 
m.eter are exposed to the same heat flux and the emf output of the test 
instrum.ent is recorded. By assuming that the standard calorimeter 
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v/as itself calibrated so that it actually indicates the true incident'flux 
and, assuming a test instrument sensor absorptivity of a = 0. 9 (Figure' 
3-22), a relation similar to the calibration Ccca be developed. ' Using a 
\'ie\v factor of unity and multiplying Equation’ 3-36 by the absorptivity 
gives the laeat flu>: absorbed by the instrum.ent sensor during calibra- 
tion as • ■ 

>\= 5.6 

r 

n 


- Sen ~ ^S ' eff 


(X) d\ 


(3-3S) 


\= 0.2 


and 


= 0.S4, crs = 0.9, «lA, -Sen =. 


X=5. 0 


Ei.(X) dX 


X=0,2 


v/hare qA_, _Ser ^-pp^oxinaate heat'fluo; absorbed by the sensor 




during calibration. 


A-ssuming an emissivity of £ 3 ^ =1.0 for the. calibration heat 
source, Equation 3.-28 gives the incident flux indicated by the standard 
instrumient. From the latter and Equation 3-38 the calibration relation, 
which is shown in Figure 3-27, can be established. By considering 
the flu:-; absorbed by the ' sensor, both during calibration and during 
exposure in the unknown environment, the error, introduced can be 
identified and approximately corrected. 
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which gives, for the conditions 

Tri= 0. 9, Tn = 0. 84, and Tgff = 


(0 = 9oy 

pX=5. 0 

lJ 

= 0 » 9 > g ~ 0 . 9 , 

^A, -Sen = 0*711 \ Ej (X) dX 
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During calibration the test instrument is exposed to a heat source 
v/iDi a view factor of unity. Using the standard instrument's indication of. 
nic incident flm: and the measured ernf,. the ilmc absorbed by the sensor 
is related to that incident on the Avindov/* During exposure to an uihmown 


flux, the process is reversed. 


The fiucc absorbed by the instrument 


gcnersrccs an emf which is then related to the incident flux indicated 
duri::g calibratipii. An error can be introduced because the unlmown 
incident flux is equal to the inc.ident flux during calibration only when 
the conditions are identical, i. e., when the emitted flux is of the grey- 
body type and the view factor is unity, 


For the conditions which existed during calibration, the effective 
transmissivity of the window has been shown to vary between 0. 734 ~ 

(6 - 90) i 0. 7895 .because of the angular .dependence of the window 
transmis'si\n.ty; During exposure in an uriknown environment it must be 
assumicd that the instrument will be exposed to a source which does hot 
fill the complete field of .vfew of the instrument. For this condition the 
effective transimissivity is different from that which existed during cali- 
bration and, as a result, more of the energy incident upon the window will 
arrive at the sensor. In order to correct for this difference, it is • 
necessary to know the emf- incident fime characteristics of the instru- 
ment for effective transmissivities other than that associated with the 

t 

original calibration (i, e, , 0.734 ^ ^eff ~ " 0.7S95). An approxi- 

mate correction .procedur.e can be developed by determining how the 
incident flux. would’ have varied with changes in the effective transmissivity 
while holding the flux absorbed by the sensor constant. 

3y integrating Equation 3-37 over limits which range from 
6, = 0 to 92 , the effective transmissivity as a function of ^ can be 
obtained. Thus 
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Rewriting Equation 3-4i gives 


p\=5. O '. 

\ (\) d.\ 

^>.= 0.2 


(q 


R, -SerJ 'ca] 


test 


"^eff (4=2 ) 


v/ner e 


=5. 0 


E, (X) dX 


(3-42) 


is the Eux over the interval 0. 2 = X 

= 0.'2 -=test 

= 5. 0 that would be incident upon the window of an equivalent ins.trument 

v/ith reff(q 2 ) during calibration to produce the same absorbed flux. 

Knowing the incident flTix over the interval 0.2 S X = 5. 0, it is then 

possible to compute the incident fl'ox over all wavelengths using Planck 

functions, tiius ■ • ■ 


qn 


.X=oo 


X=0 


Ei(X) dX 


test 


.X=0.2 


X=0 


El (X) dX 


X= 5 . 0 


El (X) dX. 


X=0.2 


A=oo 

+ \ El (X) dX 

\= 5.0 


(3-43) 


test 


To assist in evaluating the incident flux from Equation 3-43, the Planck 
function over the intervals 0 = \ ~ co and 0, 2 = X ^ 5 is plotted in 
Figure 3-29. , • 
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A sample problem can best illustrate how the indicated flux can 
be corrected* Support, for the disk of Figure 3-30, with o 2 =-25 
degrees, the.emf output of the instrument indicated an incident flux of 
iO 3t'a/ft’ - sec* From Equation 3-38 or Figure 3-27 the true absorbed 
flm: for the calibration condition would have been 


' n 


0. 84, 


-Sen) cal = ^ Btu/ft^ -sec 


and 


^n = 0-9. , (qA,-Sen)cal = 5. 2 Btu/ft^ -s ec 

By assuming that the emf output of the instrument is proportional to 
the absorbed flux, the flux, absorbed by the sensor, whenc^j = 25 degrees, 
is the same as the above. 

For 92 '= 25 degrees, the effective transmissivity of the window 
obtained from Figuf-e 3-28 is , ■ • 

Tj,' = 0.84,, . -Tgff (92 ='25») = 0.84 

or , ' • ' 

= 0.9 . ■ (92 = 25») *."o .9 . ■ . ' 

Then, using Equation 3-42, the true' incident radiation function gives 



0.84, j El (X) dX 


L^\=0.2 


Jtest 


4. 8 

0. 9 X 0. 84 


■= 6.35 Btu/ft^-sec 
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= 0.9„, 


X=5. 0 


X=0.,2 


'-'i (\) d\ 


test 


0. 9 X 0. 9 


=, 6.42 3tu/ft^-sec 

/ 1 3 . 0 

From the plot of. the radiation function, \ (X) d\ (shown 

‘^X=0.2 

in Figure 3-29), the correct incident flux can be obtained: 

Tn = 0.84, » 8.9 Btu/ft^-sec 


T^, = 0. 9 , ,qir « 9. 0 Btu/ft^-sec . 

In this example, the error in the indicated incident flux is 

= 0.84, • e = 100 X — ■ g" = 12.4% 

= 0.9 , e = 100 X =11.1%,. 

Following the same procedure as indicated in the preceding 
example, the ratio of indicated incident flux to actual incident flux for 
various included half angles was computed. (The results are plotted 
in Figure 3-31. ) As seen in the figure, the difference 'between the 
indicated flux and the actual flux is highest when the included half is ■ 
less than 25 degrees and increases slightly with the incident flux. 


The included half angles of the Saturn S-IC heat shield radiation 
calorimeters were evaluated and are shown schematically in Figure 
3-32. The maximum included half angle of the heat shield radiation ' 
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INDICATL'D INCIDENT FLUX 















FIGURE 3-32. S-IC HEAT SHIELD RADIATION CALORIMETER INCLUDED HALF ANGLES 




caiorirrieter s was approxir -ited to be 30.5 degrees at sea level. Because 
^of the orientation of the calorimeters with respect to the engines, signi- 
ficant plume omansion would be required to increase the included angle 
beyond 30. 5. degrees; therefore, it must be concluded that tl-ie S-IC 
indicated radiation heating rates probably were in error by at least 11 
to 12 percent throughout most of the flight. 

An additional error could have been produced by the nongrey 
plume radiation. Figure 3-29 shows the blackbody radiation contained 
within the wavelength range of 0. 2 micrometer = X = 5. 0 micrometers 
compared to the corresponding blackbody incident- flux. If the plume 
radiation was not continuous but emitted radietion with significant 
spectral windows characteristic of radiating gases, then a significant 
error could result. Considering the previous .example, .when the inci- 
dent flux was 10 Btu/ft^-sec, the radiant energy contained within the- 
wavelength range 0. 2 micrometer - \ = 5.0 micrometers was 6. 35 
Btu/ft“~ sec (with Tyi - 0. 84). If the disk w^ere such that it radiated 
energy in the amount of 6.35 Btu/ft“~sec only within the range 0. 2 
'microimeter = X = 5.0 micrometers, the radiation calorimeter would 
respond approximiately the same as if the disk were a blackbody radiat- 
ing at the rate of 10 Btu/ft^-sec. As seen by the difference, this could 
product a significant error in the indicated results. 


3.5- BASS PRESSURE PROBES 

3. 5. 1 Physical Descriptions 

The Saturn base pressures probes consisted of a tube connected 
to an electromechanical transducer. The transducer translated the pres- 
sure being sensed by the tube into an electrical signal. All pressure 
instrutTients were calibrated prior to flight use. 
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3.5.2 Base Pressure' Ivleasur ernent Errors 

The measured pressures are another source. of. error* When 
the raw flight data 'is reviewed, it is. seen that the pressure probes 
indicate pressures prior to vehicle lift-off and after engine cutoff which 
do not agree with the known ambient conditions. In some cases, negative 
absolute pressures are indicated after engine cutoff. These measure- 
ments are well within the specified accuracy of the instrumentation, but 
are not suitable for determining the base flow Mach numbers without 
correction. 

Figure 3-33 shows a typical S-IC base pressure measurement. 

At the top of the figure is shown the raw measured data obtained from, 
measurement D36-106, AS-502 flight. As seen in the figure, the 
indicated pressure is too high prior to lift-off and too low after engine 
cutoff. By averaging the naeasured. values taken from the digitized, 
printout of this measurement, it can be shown that the indicated value 
prior to engine ignition was approximately 15. 025 psia, whereas the 
ambient pressure was 14.75 psia. Likewise, the indicated value at 
130 seconds was approximately -0M30S.psia, whereas the ambient 
pressure at this time was 0. 0276 psia. The specified accuracy of this 
instrument was ±3 percent of the full s cale range which gives a possible 
•error of ±0. 6 psia. By comparing the indicated error .(i. e. , 0. 275 psia • 
at lift-off) to the specified accuracy, it is seen that the actual error is 
‘ much less and, although this is encouraging, does not eliminate the 
probability of .analyzing the flight data. 

One of the basic parameters of interest in analyzing the. base 
convective heating is the base flov/ Mach number. In subsonic com-, 
pressible flow, the Mach number can be computed if the static and 
total pressure (T^.S> specific h<^at ratio [y) are known. 

The equation which gives the Mach number is 
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PRESSURE (psia) 




imiEHT PRESSURE AT LiFT-OFF = 14.75 psia 
MEASURED PRESSURE PRIOR TO LIFT-OFF (D36-.106) = 15.025- psia 
DELTA PRESSURE = 0.275 psia 
. MEASURED PRESSURE AT 130 sec = - 0.1308 psia 
AMBIENT PRESSURE 130 sec = 0.0276 psia' 
lECO = 144.39 sec 
OECO = 147.34 sec 

AMBIENT PRESSURE AFTER ENGINE CUTOFF = 0.006 > P > 0.003 psia 
FIGURE 3-33. AS-502 BASE PRESSURE MEASURED AT LIFT-OFF AND CUTOFF 
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In base flow the specific heat ratio is knowca since the gas is primarily 
ambient air (y = 1*4) during the aspiration phase and engine exhaust 
products (v 1.2 for LOX-RP-i) during the flow reversal phase. 
Therefore, if the total and static pressures can be measured during 
the flight, the Mach number can be evaluated, approximately, iising 
Equation 3-44. ' For this purpose, static and pitot probes were located 
in the base region-of the S-IC vehicles. The pressure measured by 
one set of pitot and static probes during the AS-501 flight is shown in 
Figure 3-34. The ambient pressure and the experimer.tally obtained 
pressure in the base of a bluff body are also shown in the figure. If 
these measured values are used in Equation -3-44, they indicate a 
base flov/ Mach number of 0. 1 to 0. 2 shortly- after ‘lift-off; however, 
the indicated f low direction is out of the base, which is opposite to the 
un-derstood characteristics of the base flow. This trend also persists 
for a significant time over the first phase of the flight. Obviously, 
the raw measured pressures are not accurate enough to allow the 
computation of the Mach number and flow direction; this situation 
can be improved, however, if the measured values can be corrected. 
Since the actual pressure is knovm (approximately) prior to engine 
ignition and after engine cutoff, these two end points can provide a means 
of correcting the 'measured values. The correction technique is shown 
schc'matically in Figure 3-35. The technique consists of determining 
the difference between the measured value and the ambient pressure 
prior to lift-off and after engine cutoff. If the indicated error is plotted 
versus the measured pressure and a linear relation is assumed between 
the end points, the measured values can be corrected. The measured 


3-74 


m 


i I 


n 

LJ ' 


a! 


i i : 


I 

rn 

1 ( 

l_ - i 


U [ 


u n 


ri ! 

U i 


i i f 


u 


r J ' 

U 1 


n i 


f j i 


n ' 


i r 




MEASURED PRESSURE CORRECTIONS SCHEMATIC 



pressures shown in' Figure 3-34 were corrected in this way and are 
showii in Figure 3-36, Because of the nature of the tecluhque, the final 
results are only approximately correct; however, the obvious errors 
have been eliminated.’’ 

The base flow Mach number can be approximated from tliese 
corrected pressures. by assuming that the maximum absolute pressure 
is the total pressure and the minimum pressure is the static pressure. 
The computed results are shown in Figure 3-37. The Mach number, 
computed using the bluff-body base and ambient pressure, is also shown 
in this figure. ■ ' ' 


3.6 CONCLUSIONS AND RECOMMENDATIONS REGARDING 

INSTRUMENTATION 

The results of this study have produced the following conclusions: 

o All instruments used to measure the Saturn inflight 
base heat transfer parameters ■ appear to have had a 
significant potential for erroneous indications.' 

o The gas temperature probes were sensitive to. radia- 
tion heating and lacked response in a low-mass flow 
rate environment; as .a result,- the indicated gas 
temperatures could be higher or lower than the actual 
gas temperatures. 

o The total heat flux meters were sensitive to the mis- 
. match in surface temperature between the sensor and 
the surroundings ; this could have- produced a high indi- 
cated convective flux. Condensation of mpisture on 
the surface of the Instrument could have produced a low 
, , • indicated convective flux. ■ ^ . ■ 

G The radiation calorimeter window combined v/ith the 
^ calibration results to produce indicated incident heat- 
ing rates which were higher than the actual flux. 

o The measured . results indicate the masfnitude of the 
individual base heating parameters but, to an unknown 
and indeterminable accuracy. 
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FIGURF 3-36. 'AS-501 CORRECTED BASE PRESSURES. 





501 i-LfGHT RESULTS 

502 FLIGHT RESULTS 

USING BLUFF BODY AND 
AMBIENT PRESSURE 
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coir/bined without liberal modification, to delineate 
important aspects of the’base heating problem 


o The results point up a. serious need for the develop- 
ment of accurate gas tenrperature probes and heat 
flux meters which can be used in an environment 
similar to that of the Saturnbase region, 


V/ith respect to future flight test instrumentation, the followmg 
commendations are made: ■ 

o For booster stages where the radia.nt flux is signifi- 
cant and the mass flow flux probably changes by 
several orders of magnitude during the flight, it is 
doubtful if an accurate gas temperature probe can 
be developed -using a conventional thermocouple. 

However, inore accuracy and greater confidence 
in the results can be obtained by calibrating the 
instrument in a combined radiant and convective 
environment and using the calibration results for 
evaluating the approxi-mate base gas recovery tem- 
• perature. For upper stages where the radiant flux 
and mass flovy flu:<’are both low,, the probe should. - 
also be calibrated because the*, reradiation from, the ■ 
probe junction may still cause- erroneous indications. 

o For the total heat flux meters/ a properly designed ring 
which surrounds the instrument at the exposed surface, 
could reduce the' effects of moisture condensation and 
surface temperature mismatch. This instrument should 
also be tested in a combined, radiant and convective ' . *. . 
environment to de.termine if empirical correction pro- 
cedures are necessary. 

o For the radiant flux meter it is doubtful if an instrument 
can be developed which does not have errors caused by 
variations in the window transmissivity. A method of 
approximately correcting the measured results has 
been outlined in Section 3.4. This method should be 
verified experimentally and on future flight vehicles the 
instruments should be oriented with respect to the 
radiant source so that a correction procedure can be 
used* 
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i or pressure measuretTients, the existing instrumenta- 
tion IS probably as aecurate as Ccor reasonably be 
e:-rpecteci. Therefore, for future flight test's 'it is 
recorrnnended that the same type instruments be used 
and corrected as outlined in Section 3.5. . ■ 


} i 


n : 




n il 


■n ’ 
( ! ; 


u 


P7' > 

i. ; 


Li i 


LJ : 


< — j 

i 

LJ ^ 


f ^ I 

t i i f 



REF 


:RENCES - SECTION 3 


Eencdict, B. U., et ai, "Eval-aation Test of Thermocouple,, 
Tyoes 9 and .35, OR.D P/Ns S9600S5 and 50M10099"> Chrysler 
Corporation Technical Memorandunn EER-147, Prepared- 
under Task Order GS;C-C-29, Contract NASS-6, December 
1961 


3-2. Connell, H. A., "Evaluation of Saturn SA-1 Base Heating", 
NASA Report MTP-P&cVE-P-6 1 -21, December 1961 

3-3. "Qualification Test Report for Sensor (Transducer) Gas 

Temperature", Boeing Document REC46727A, Rev. A, Attach- 
ment to 5- 1 1 1 8-M-7 1-342, August 1967 (Qualification Test 
• Report on S-IC Engine Gas Temperature Probes, retained in 
MSEC Document Repository) 




3-4. "Report of Test on Transducer, Temperature, .Ambient Gas", 
Boeing Documient QTR-33F, Revision A, June 1966 (Qualifica- 
tion Test Report on S-IC Heat Shield Gas Temperature Probes, 
retained in MSFC Document Repository) 

3-5. ■ Sudunas, j. A., "AS-503 S-II Base' Region Gas Recovery 

Temperature Evaluation", .North American Rockwell Corpora- 
tion, Internal Letter S-II- 1 9P-700-69-2 1, March 1969 

3-6. Sudunas, J. A., "S-II Base Region Gas Recovery Tem.perature 
Probe Analysis Programij YF.0028", North American Rockwell, 
Internal Letter LEVA-190-405-70-35, June- 1-970 

3-7. "Design and Development, Testing of Transducer Temperature 
Thermiocouples V, RdF Corporation Developm.ent Test Report No. 

. 91, April 1966. (Test -Report on S-II Stage Gas Temiperature 
- Probe) - - ' . . 

3-8. "Major Design Review Supporting Data, .North Am.erican 

Aviation, Inc., S&ID Procurem.ent Specification MC449-00.06, 
Thermocouple Temperature Transducer", RdF Corporation, 

May 31,' 19*66 ■ ' * . ' 

3-9. Wiebelt, J. A., Engineerin,g Radiation Heat Transfer , Holt, 
Rinehart and Winston, Inc., 1965 


3-81 


r 

i' 



{■ 

.1 


! 

\ 


!. 


i 


I 

1 

i 

I 

» 


.1 




I 1 ■ 

\ I 



f 


I 




Ulo. 


REFERENCES - SECTION 3 - Concluded 


McAdams, V/. K. , I-Ieat Tzmnsmission , 3rd Edition, McGraw- 
Hill Publishing Company, 1954 


3-1 i. Buchrnann, R.. C., J. T, Cbamibers, and W, H. Giedt, ’'Invcsti- 


f/ation of Surface H 


Tux Me 


with Calorimeters*', 


3-12. 


ISA Transactions, pages 143 through 151, April 1965 

Bender, R. L. and W. B* Shaffer, "Examination of Saturn V 
' First Stage Thermal En\dronrnent Flight Measurements 501- 
511", Boeing Company Document 5-9410-H-44S, September 
1972 


3-13. 


1 A I ' 


3 - 15 . 


3-L 


Thermal Instrumentation", Hy Cal Bulletin, No, A-102, 1970 

"Properties of Infrared Transmitting Material' 

Bulletin, Servo Corporation of America, 1961 

Garbuny, M, , Optical Physics , Academic Press, Inc*, 1965 


Li 

n 

LJ 

n 


c > 

u 


ecliniques", 

Pi 

Li 

March I960 

A-102, 1970 

t — ^ 

u 

Servo /Infrared 

n 

LJ 

Inc., 1965 

f 

■ ' - i 1 


\ \ 
f ) 

\ 




1 


n 


[ 1 


! — ! 





3-82 






4. FLIGHT TEST RESULTS ■; 

4. 1 BASE :-:>:YIRONME1nT measured during SATURN 

FLIGHT TESTS . , ‘ • 

The Saturn flight test data have been reviewed with regard to 
the evaluation of the in-flight base. heating. . As anticipated, a volumi 
nous amount of Saturn flight test data was available for evaluation. 
The measured data was obtained from the following sources: 

o Raw Flight Data - This data, consisting of plotted or 
tabulated time histories of the measured parameters, 
was obtained from the MSFC Computation Laboratory- 
flight data microfilm laboratory.’ '■ In general, this 
data was used in the study where possible. 

o Contractor Flight Reports - For some of the flight, 
tests, the stage contractors published separate aero- 
thermodynamic flight evaluation reports. Some of 
these data were used in the evaluation, primarily, 
where the raw flight data was not available. 

o Original Flight Evaluation Working Papers - These 
data consisted of raw flight data, as well as calibra- 
tion curves and original flight evaluation not pub- 
lished in other reports. These data were obtained 
from the MSF C Astronautics Laboratory, Thermal 
Engineering Branch files. 

o FEV/G Reports - These are flight test reports pub- 
lished by the MSEC Flight Evaluation Working Group. 
These reports contain a summary of the vehicle's 
flight performance and measured parameter s. 

o .Trajectory Reports - These reports contain the tra- 
jectory of the flight vehicle and ambient conditions at 
or near the time of launch. These reports are main- 
tained in the MSFC Aero- Astrodynamic s Laboratory 
files. 
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o I? and C Lis* - These are documents which list the 
tyne anci location ci ail instruments that are on trie 
vehicle, r/ost of these documents were available 
in the I'/SFC Document and Drawing Repository,. 

A review of the flight test d^ta indicated that it would be beneficial to 
plot the measured results of each mieasurement for all flight tests- on 
a single plot. This v/as found to be beneficial because any given param- 
eter rr:easured at the same location in the base region was found to 
vary siighAiy from flight to flight. For some flight tests.these varia- 
tions could be att .buted to minor variations in vehicle configuration. 

For most, however, no obvious reason for these random variations 
appeared to exist.-. Therefore, plotting the data measured at one loca- 
tion during several test flights -could be used to identify measurements 
which probably failed and/or those which indicated excessively high or 
low values during the flight. 

For the upper stages the measured data could be -plotted directly 
from the rav/ flight test results. For the booster ■ stages, however, it 
was necessary to normalize the results. This;. was found to be necessary 
because most flight test vehicles flev/ slightly different trajectories. 
Because the vehicles were at different altitudes., at the same flight 
time and because a miajority of the base heating parameters measured 
were primarily a function of local ambient pressure, these differences 
had to be accounted for to obtain consistent data for comparison. As a 
means of normalizing this data, all of the measured results were related 
zo the Saturn V, AS -501 trajectory by using altitude as the common 
factor. In this way, all of the booster flight test data could be plotted 
as if the vehicles had all flown the A.S-501 trajectory (for more details 
of the normalizing procedure see Reference 4-1). This was found to 
be a painstaking procedure, but resulted in data that could be comipared. 
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All of 
above 
da^a I 


the imoor tant flight test dat 
and is presented in Volume 
or each stage is presented. 


:a was. plotted in the manner described 
II.' r’or this report, only representative 
For an alternate analysis,’ see 


• * * Reference ,4-2. 




4.2 ANALYSIS' AND COMPARISONS OF SOME OF THE MEASURED ' 

SATURN FLIGHT TEST DATA 

4.2. 1 Saturn I Flight Results 

Figure 4-1 shows the bands of gas temperature measured in . • 
the base region of the Saturn I, block II vehicle.’ • The apparent times of 
the start of flew reversal are also shown in this figure. The time of 
■ flow reversal can only be approximately established and the particular 
times shown in Figure 4-i were .chosen because of significant increases 
in the measured temperatures. 

Typical pressures nieasured in the base region of the Saturn I, 
.block II vehicle are sho'-n in Figure 4-2. In this figure, it can be seen 
that the heat shield pressure is very nearly the same as the local ambient 
pressure. It is also seen in this figure that the flame shield pressure 
is significantly different from the ambient or base pressure. Since 
pressure is indicative of the flow field, these measured pressures 
appear to indicate. a weak flow field in the region of the heat shield and 
a strong flow field in the flame shield region. Another factor which is 
indicated by the flame shield measured pressure is the probability of a 
so-called "choked*' flow condition (i. e. , condition where the mass flow 
rate is unaffected by any further decrease in the ambient pressure). 

When rhe choked condition is reached the temperature (see Figure 4-1), 
pressure (see Figure 4-2), and heating rate (s.ee Figure 4--3) become' 
constant. 
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Typical £larr:e shield heating rates measured on one of the 
Saturn I,' block II flight tests are shown in Figure 4-3. The total heat 
flu:*: includes both the radiant flu:*: (also shown in the figure) and the 
convective flu:*:. The difference between the total and radiation flux 
gives the convective flux. The results shown in Figi-re 4-3 indicate 
that a significant amount of convective heating occurred both prior to 
and after the apparent start of flow reversal. -By observing the results 
of this figure, it can be seen that the convective- heating (i. e. , the 
difference between the' total and radiant flux) was nnuch higher after the 
flow reversal than prior to flow reversal. However, the fact that the 
measured gas temperature is high and the convective heating starts at 
lift-off and continues over the first half of the flight suggests that sooie 
of the high temperature exhaust gases are being reversed even at lift-off. 
Figure 4-3 also shows that the radiant flux is the only mode of heat 
transfer after. 60 seconds. This might not be a completely accurate 
conclusion since there were problems encountered in evaluating the 
total heat flux which are explained in Reference 4-3. However, it does 
appear that radiation was the dominant mode of heat transfer. 

Another interesting characteristic can be observed in Figures 
4-1 and 4-2. During one of the flight tests, a center engine was cut off 
appro>:irriately 20 seconds earlier than the other three engines. .The 
resulting reduction in the flame shield temperature and pressure are 
shown in the figures. The gas temperature dropped approximately SOO^^F 
and the pressure dropped approximately 2 psi. The heat flux (not shown) 
dropped only slightly and the choked flow condition appears to have been 
r eestabiishec: . ■ *. ■ 

Shown in Figures 4-4 and 4-5 are representative bands of 
radiation and total heat flu>: measured in the outer and inner zones of 
the Saturn I heat shield. In the figures, it can obviously be seen that 
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FIGURE 4-5e SATURN I, BLOCK II INNER REGION HEAT FLUX 




radiation was the cornina*nt source of hear input Throughout the flight. • 

By comparing the total heat flu>; and the radiant flux, it can be seen 
that a significant amount of convective cooling exists, -particularly during 
the early phase of the flight. 

4.2.2 Saturn S-IC Flight Results 


One of the main advantages of normalizing the flight test results 
is that it allows direct comipariscn of the same data for different vehicles. 


Figures 4-5, 4-6, and 4-7 show comparisons between Saturn I and Saturn V, 
S-IC base therimal environment. Figure 4-6 is a comparison between the 
Saturn I (outer region) and the Saturn Y, S-IC base gas temperature. 
Surprisingly this shov/s the gas temperature and the apparent time of 
flow reversal to be similar even though significant configuration dif- 
ferences existed between the two vehicles. 


Figure 4-7 shows a comparison of the radiant flux measured in 
the base region of these two vehicles. The radiant flux incident upon the 
3-IC heat shield was much lower than that incident upon the Saturn J 
heat shield. This difference was primarily due to the relative location 
of the heat shields and the location of the engines with respect to each 
other. The heat shield on the Saturn I vehicle W'as located approximately 

2. 4 engine radii above 'the exit plane, whereas the S-IC heat shield was 

3. 26 engine radii. The Saturn I engines were spaced farther apart, 
with respect to a point in the outer region of the heat shield, than the 
S-IC engines. Both of these factors tended to increase the radiation 
view factor (see Equation 2-17) between the plumes and the heat shield 
resulting in higher incident radiation to the Saturn 1 heat shield. 

Another interesting characteristic of the S-IC radiation is the 
sharp increase ’measured after the apparent time of flow reversal. 

This same characteristic was evident in the Saturn I flame shield 
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region (see Figure 4-3) and, to a lesser degree, in the heat shield 
region. This is thought to be caused by the start of full flow reversal 
and will be discussed later. 

. ■ A comparison of the Saturn I and S-IC total heat, flux is shown in 
Figure 4-8. The uotal heat fTux is somewhat indicative of the overall 
thermal environment since the area under the curve represents the net . 
energy input to the surface. In light of the latter, a comparison of the 
Saturn I and S-IC total flux of Figure 4-8 would indicate that the S-IC 
heat'shieid environment v/as slightly less severe than the Saturn I. 

In convective heat transfer, the flow velocity is an important 
parameter. In compressible flow the velocity can .be computed if the 
.gas temperature, total pressure, and static pressure are known (see 
Fquaticn 2-9)- In order to measure the static and total pressure, pitot 
static tubes were located in the base region of the Saturn S-IC vehicle. 

The results measured at t\^ ^ locations in the S-IC base region are shown 
in Figures 4-9 and 4-10. The measured results . shown' in the figures 
indicate that the freestrearn ambient air would flow into the base region 
during the first half of the flight (i. e. , the bas e 'pres sure being lower 
than t?:e ambient' pressure indicates inward flow). Likewise the base 
pressure being greater than the ambient pressure after 75 seconds indi- 
cates reversed flew back into arid out of the base region. The evaluation 
of the base flow velocity will be discussed later. 

4. 2. 3 Saturn S-IV Flight Test Results 

Tv/o of 'i:he Saturn upper stages had clustered engine arrangements 
and both used LOX/hydrogen as the propellants.. The Saturn I upper stage 
consisted of six RL-10 engines clustered’ about a common center. The 
instrumienfation in the base region consisted of black and gold-coated, 
heat flux meters and static pressure probes. Because of the LOX/ 
hydrogen propellants, the radiant flux, in the . base region was negligible • 
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■FIGURE 4-10. SATURN V, S-IC STAGE ENGINE PITOT PRESSURE 





could not be detected by combining the two different type heat flux 
meters. Engine ignition occurred at such a high altitude that the choked 
flow condition was established almost immediately. -The heat shield 
heat flux measured at various radial distances from the center are’ 


shown in Figure 4-11. Because of the very low density of the base gas 
at these extremie altitudes, the base heating rates are less than 2 3tvi/ 
ft^- sec. Another factor v/hic'O also affected the base environment on this 
stage was the' location of a heliumi heater nozzle in the center of the heat 
shield. The purpose of the helium heater was to warm the helium which 


was stored in the liquid hydrogen tank for use in pressurizing the liquid 
oxvgen tank. The heater consisted of a low temoerature combustion 


chanaber /heat exchanger combination attached to a small nozzle. The 

nozzle dis charged gas into the base region at a flow rate of 0.05 Ib/sec 

and a temperature of 800^ F. The secondary flow from this helium heater 

nozzle had some .'effect upon the base convective heating. The effect would 

appear to be a reduction in the heating rates. . ‘ • 

0 ■ 

Shown in Figure 4-12 are heat shield static pressures measured 
at various radial locations during one of the S-I'V stage flight tests. 

As shov/n in the figure, the pressures are low and indicate that the 
choked flow condition is reached shortly after ignition. 

4. 2. 4 Saturn S-II Flight Results '- 

Shown in Figures 4-11 through 4-19 are typical gas temperatures, 
heating rates, -and pressures measured. in the base region of the Saturn V, 
S-II upper stage. This vehicle had five LOX /hydrogen fueled J-2 
engines. The environment measured in the base' region of this stage 
shown in Figure 4-11 consisted of , approximately 3 Btu/ft^-sec total 
heat flux, a I 3tu/ft“- sec. radiant flux, and a lOGO^'F gas temperature. 

The resulting convective flux is, shown in Figures 4-14 and 4-15. A 
cylindrical interstage with a diameter the same as the basic vehicle 
extended approximately 7 inches below the engine exit plane. This 
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inrerstac^e was iettiscned after the first 30 seconds of the flic;ht. The 


effects of the interstage on the base thermal environment can be seen 
in the figures. The heating rates and pressures above the hear shield 
tended to be higher with the interstage on the vehicle and drooped 
drastically when it was jettisoned (see Figures 4-16 and 4-19h The 
pressure on the heat shield was also higher during this' period (see 
Figure 4-18). ' The heat shield gas temperature and total heat flux Were 
slightly lower with the interstage on than with it off (see Figure 4-15). 

Another characteristic which can be seen in -these figures is 
the decrease in heating rates and pressures associated with the shift 
in nnixture ratio.- The J-2 engine mixture ratio was shifted from the 
initial O/F = 5, 5 to an OF = 4.7 late in each flight. The decreases 
can be seen in the. figures at flight times corresponding to 400/450 
seconds. 

4. 3 ANALYSIS AND CORRELATION OF THE MEASURED 

RADIANT FLUX ' ‘ 


The radiation measured in the base region during the Saturn 
booster stage test flights w^as correlated. The correlation consisted 
of determining the ratio of rad .ation measured at any time ‘during the 
flight to the measured ground level value. The ground level value was 
chosen after the vehicle had cleared the launch pad. so that the effects 
of the flame deflector could be eliminated. The resulting ratios -are 
presented as a function of vehicle altitude in Figures 4-20 through 4-23-, 
Figure 4-20 shows the dimensionless radiation ratio obtained from the 
Saturn I measured data for tw'o different zones of the heat shield. The 
Saturn I, Block I (SA-1 through SA-4) and Saturn I, Block II (SA-5 
through SA-10) are both shown in this figure. The radiation ratio 
obtained from the Saturn S-IB flight test data is shown in Figure 4-2L 
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FIGURE 4-23. FRACTION OF GROUND LEVEL RADIATION, SATU 
AMD OTHER VEHICLES 
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Also shown in this figure is the Saturn I, Blocks I and II data band. 


studying the data, 
from, the Saturn I, 


a siigrit upward trend' appears to be indicated in going 
Block I vehicle radiation decay to the Saturn IB results. 


This trend was probably caused by slight changes in the vehicle configu- 
ration and the uprating of the engine from 165, 000-pound thrust, Block I. 
(578-psia chamber pressure) to 200, 000-po.und thrust, S-IB, (639-psia 
chaimber pressure). 


Figure 4-22 shows. the radiation ratio for the Saturn. S-I and S-IB 
flame slaieid regions. The heat shield data bands are also shown in this 
igure. As seen in this figure, the Saturn S-I and S-IB flame shield ■ 
radiation ratios are approximately the same up to an altitude of approxi- 
mately 15, 000 feet and decrease at a much faster rate than the heat 
shield radiation. Beyond 15, 000 feet a completely different trend is 
indicated. This difference can be explained and resulted from changes 
in the turbine exhaust gas discharge. The Saturn I turbine exhaust gas 
was discharged overboard through ducts on the sides of the vehicle. 

This was changed on the third Saturn S-IB and subsequent vehicles. 

For these vehicles, the turbine exhaust was discharged into the flame 
shield area. The sharp increase in radiation shown for the S-I flame 
shield coincides with the beginning of flow reversal. This same trend 
(i.e. , slight hump in the curve), but less significant, can be seen in the 
hecvt shield radiation as well. This indicates that the radiation is being 
influenced by the flow reversal. 

Shown in Figure 4-2 3 are the radiation ratios for the Saturn S-TC 
heat shield, the ATLAS, the JUPITER, the TITAN IIPC, and the Saturn - 
S-I and S-IB heat shield data bands. The S-IC vehicle had five LOX/RP-1, 
F-i engines. The ATLAS had three engines in-line, and the two outboard 




engines were the saline as the Saturn S-I engines. .The JUPITER vehicle 
'had only one engine, which was also the same'as the Saturn S-I engines. 

The TITAN III-C had two solid propellant engines. Two of the interest- 
ing trends shown in this figure are the JUPITER and Saturn S-IC results. 
Because -the JUPITER had only one engine, -flow reversal between engine 
'exhaust jets did not exist. The JUPITER radiation ratio decreases con- 
stantly throughout .most of the flight. Because the engines were'the 
same, comparing the JUPITER results to the Saturn S-i results giyes 
some indication of the engine clustering effects. The Saturn S-IC results 
are approximately the same as the ‘Saturn S-I up to an altitude of 35, 000 
■ feet. Beyond this altitude, the S-IC radiation ratio shows a sharp increase 
to a- value, well above the ground level radiation. This also coincides 
v/ith the beginning of flow reversal. 


Because of the very complex nature of the problem, accurate 
radiation predictions are very difficult to make. The technique used 
in predicting the base radiation for the Saturn boosters was a senni- 
enrpirica! method which depended strongly upon exper imental- r esuits . 
The radiation at sea level was computed using Equation 2 -17. The 
emissive powers and plume shapes were obtained from measurements 
and pliOtogr aphs taken during captive engine test firings. The radiation 
at altitude was predicted by extrapoiaiing from one class ox vehicles to 
the next using results similar to those shown in Figures 4-20 through 
4-23. If it had not been foi* using conservative estimates of the sea 
level radiation, this technique would have led to serious underpredic- 
tions of the radiation at higher altitudes. An example of this can be 
seen in Figure 4-23 by corrxpa ring the Saturn S-IB radiation decay to. the 


Saturn S-IC radiation. 



3 y c o n 5 : d e r in g Z qu a t i o n s 2 - 1 5 h r ou gh 2-17 th e pr o b 1 e m s \v i t b. 

•*:h:s to chnicue “can be discussed. V/hen the ratio of radiation-at-altitude 
CO radiation -at- sea level is plotted as- a function of altitude, this essen- 
tially shov/S' ho\v the radiation varies with ambient pressure. When this 
ratio is used to extrapolate the radiation from one vehicle to the next; the 
assumntion thai is. being made is that the plume emissive power and 
view factor will vary with pressure in the same general way. This, 
in ‘:urn, is' assuming that the plunae temperature, emissivity, and shape 
will vary with ambient pre s sure in this same general way. By consider- 
ing Figure 4-23, this is obviously not a valid assumption. By comparing 
the radiation character istic s of the JUPITER and Saturn S-I and S-IB 
(Figure 4-23) and the Saturn S-I flame shield (Figure 4-22) where the 
engines were nearly the same, it is seen that clustering has a significant 
effect upon the radiation at higher altitude. Furthermore) considering 
ihe Saturn S-I flame shield radiation when liow reversal starts 
(approximately 20, 000' it) and Equation 2-17, it appears that the after- 
burning is either flushed back into the flame shield or' starts to burn 
nruch closer to the engine exits at this time. In comparing the Saturn S-I 
flame shield radiation to the Saturn S-IB with the flame shield turbine 
exhaust discharge, it would appear- that dumping this relatively large 
amount of unburned fuel into this region prevented the afterburning from 
taking place. In comparing these results to the S-IC radiation shown. in 
Figure 4-23, it appears that afterburning and flow reversal effects also 
caused the increase in radiation on this vehicle., In comparing the 
Saturn S-IC radiation ratio to the Saturn S-I heat shield radiation in 
Figure 4-23, where the clustering arrangement is not too dissimilar, 
it must be concluded that afterburning and flow reversal effects take ' 
place in .very complex ways. It must also be concluded that extrapolat- 
ing the radiation from one vehicle configuration to another should be 
approached with caution. 
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4. 4 S VALUATION AND COMPARISONS OF THE SATURN S--IC 
. CONVSCTIVE HEATING 

•Cor relating a^nd normalizing experimental data is a necessary 
step in, solving a convective heating problem which has hot been previously 
solved. However, the serious instrumentation-related, problems associated 
with much of the Saturn flight data have made the correlation of the flight 
results very difficult. Such results as gas temperatures indicating con- 
vective heating simultaneously with the heat flm: meters indicating 
cooling are factors which reqhire serious explanation in correlating 
the flight data. The S-IC flight results are typical of the problems which 
are encountered. 

V 

■ Figure 4-24 shoves the base gas temperatures measured during 
a ‘typical S-IC flight test. The results shown are from probes which 
were mounted through the heat shield with the sensing element located 
from 0.25 to 2.5 inches below the heat shield surface. Also shown 
in the figure is the center temperature of the- total calorimeter, membrane. 
The mernbrane center temperature was computed from the measured 
total heat flux and the calorimeter sink temperature using the following 
equation. 


= qt 


4K 6 


(4-1) 


wnere 


foil calorimeter center temperature 
foil calorimeter sink temperature 


radius of the foil 


conductivity of the foil 
foil thickness. 
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Also shewn in the figure is the total (i. e. , stagnation temperature of 
the free stream). 

Shown -in Figure 4-2 5 are the bands of radiant and total heat ihix 
measured by the one set of heat shiehd calorimeters during the first five 
S-IC flight tests. Included In Figure 4-2 5 is the adjunction (denoted by 
plus signs) in the radiation band which would account for the error in 
measured Tcidia.tion resulting from angular variation in window trans- 
missivity (see Section 3.4,2). 

Shown in Figures 4-26 and 4-27 are the convective heating rates 
(with some of the obviously erroneous data deleted) evaluated from the 
heat flux meters at two locations on the S-IC heat shield. The convective 
flux is shown for both the uncorrecred radiation (i.e,, as measured) and 
for the measured radiation corrected for angular variation in window 
transmis sivity. 


3y comparing the results of Figure 4-24 to Figures ' 4-2 5, 4-26, 
and 4-2 7 during the aspiration phase of the flight, the disagreement in 
measured data becomes obvious. Figure-4-24 shows the. measured 
gas temperature to be higher than the highest temperature of the- 
calorimeter beyond 40 seconds. Prior -to 40 seconds, the measured gas 
temperature is approximately equal to the calorirneter ternperature. 

The implied result would be convective heating of the calorimeter at 
least beyond 40 seconds. However,- the result obtained by evaluating 
the convective flux using the measured heating rates is just the opposite, 
i.e, significant convective cooling during the first 75 seconds of flight • 
(see Figures 4-26 and 4-27). 
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i'his disparity must be explained. before it is possible to correlate 
flight test data. Some possible reasons are: 

o Because of the sensitivity of the gas temperature 
probes to incident radiation, the temperature 
measured by the probes was higher (during the 
aspriation phase) than the actual base gas temperature. 

c If the incident radiation beyond h - S^pm were negligible 
{i. e. , plume radiation significantly different from a 
greybody), then an additional error in the measured 
radiant flux would exist (see Section 3.4.Z).- 

c The possibility ox moisture vaporized from the heat 
shield condensing oh the total calorimeter surface 
exists and could have ..caused a low indicated total 
flux (see Section 3.3.4). 


The possibility of explaining the disagreement on the basis of the 
surface tenaperature misnnatch has also been suggested. This reasoning 
has been rejected. Using the semi-empirical solutions' which 'have been 
derived for computing the heat transfer coefficient across a surface 
temperature discontunity, negative heat transfer coefficients can be 
computed for specified conditions. However, negative heat transfer 
coefficients cannot exist in reality b^jcause of violating the first law of 
thermodynamics . The situation which probably exists is that the recovery 
temperature of the boundary layer fluid which influences the heat transfer' 
near the discontinuity is also changing. These changes in recovery 
temperature 'nave not been included in the solutions of the surface 
temperature mismatch heat transfer. If this could be included in the 
solutions, negative heat transfer coefficients s'hould not be possible. 

Of the reasons proposed, none can be conclusively ruled out. 
Ho'wever, available spectral radiation measurements [Refs. 4-3 and 4-4 
(also see Figure 6-13b)j indicate that large LOX/RP-i fueled rocket 
engine plumes do tend to emit greybody type radiation. Condensation of 
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'-•TiOisture- or. the surface of the calorimeter is apparently a real possi- 
bility; however, the degree to which it could have reduced the nioasured 
toiai neat iIum is soecuiative. 

T>- desired goal of any correlation of convective heat transfer 
is the development of a dimensionless Stanton number, Prandtl number- 
Pvcynoics number relation. These correlations usually provide the 
basis for selecting the coefficients C, m, and n and the characteristic 
dimension, of Equation 2-13. Another advantage of this correlation 
is that it provides a direct means of conrparing the convective heating 
for different geometric bodies and flow fields. A dimensioi^less 
correlation of the Saturn S-IC base heat shield convective heat tranfer 
during the aspriation phase of the flight v^as developed using the following 
assumptions: ' . ■ ' 

o The ireestream total temperature (see Figure 4-24) 

^vas used as the base gas recovery temperature. 

o The transport properties of air (i. e. , p, k) 

were used and evaluated at the film temperature L 
[ 1 . e . , V t ^ 


n 


T^)/2l . 


o The base flow velocity was evaluated .as discussed 
in Section 3.5.2 (see Figure 3-35). 

o Because the flow was into the base region during this 
phase of the flight, the characteristic dimension was 
taken as the distance from the edge of the base .'to the 
heat flux meter. 

o The convective heat transfer was evaluated every 5 
seconds for each heat flux meter and the results 
joined by a straight line in. Figures 4-28 and 4-29.. 

Before correcting the measured radiation, most o'£ the Saturn 
jsults appear to fit between lines A and B shown in Figure 4-28. 
The resulting correlation eauations obtained for Bines A and B are: 


0 


f ! 


I i 


n 

t I 


0 


u 


O' r'' •>. O i 


U 


Li 


n 


Li 

n 


r j 
{ I 


\ \ 













4-45 
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transi 


data rnigrit. oc justii* 
rature ntismatch. 
•pical lor membrane 
cictor of up to 1.5. 
izec heat shield moi 
'er sh6v;m for the fli 


led by the calorimeter -heat shield surface 
If the measured'' re suits shown in Figure 3-16 
calorimeter s , the flight data miight be reduced 
Condensation heat transfer caused by the 
.sture could also account for the- higher heat 
ght data. , ■ 


Because cooling exists during the aspiration phase of the flight, 
vise of rhe correlation equations for a flat plate -would give conservative 
heat sliield design data. 








CONCI,U-SIQN-S REGARDING THE FLIGHT TEST RESULTS ' 
Some of the conclusions regarding the flight test results are as 




o Useful qualitative trends and characteristics can be 
obtained from the flight test data by comparing simi- 
rar paraineters over a series of flights and between 
similar vehicles. 




o The attainment of quantitative values and comparisons 
directly from the measured flight data should be con- 
sidered with caution in light of the instrumentation- 
related errors and problems. 

o Because all necessary convective heat transfer parame- 
ters must be evaluated and included and because the 
results can be compared directly to similar convective 
heat transfer problems, the development of the Stanton-r 
Prandtl-Reynolds number relation should be the ultimate 
goal of any base-heat transfer evaluation. At present, 
the inability to measure the necessary parameters 
a.ccurateiy on flight vehicles makes this correlation 
difficult and requires liberal interpr etation of the 
measured data. However, until these correlations can 
be made using flight-related data, serious questions 
regarding the measured data will continue -to- exist. 
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The rnea.sured radiation results for vehicle s using 
. hydrocarbon and solid propellants appear to indicate 
the existence of an empirical relationship between 
the reduction in radiation with decreasing ambient 
pressure. This pattern may offer a simple inethod 
for making good engineering estimates of the inflight • 
radiation from measurable ground level conditions. 

o Since the most serious 'problems with the flight test 
results appear to be those associated with the instru- 
mentation errors, additional effort could be well 
spent in this area. In light of the Space Shuttle vehicle, 
the iTiost pressing need is for the development of 
accurate, reliable base parameter measurement-, 
techivlques. The alternatives are 

^ Spend the time and imoney necessary to develop 
new instruments with inrproved accuracy 

A Calibrate (i. e. , to the degree possible) the exist- 
ing instrurnentation and evaluate the base thermal 
parameters approximately. 


5. SCALE /ViODEL TESTING OF BASE HEAT 'TRANSFER 


D. 1 


INTRODUCTION 


In addition to the problem areas associated with rocket exh5.x:st 
plume /vehicle inte ractions is the requirement for knowledge of the 
complete integrated system performance characteristics throughout 


the flight trajectory. Such information is often required both during 
the desig stages for purposes of optimizing new vehicle designs and 
later during flight test phases to account for variations between observed 
a.nd predicted performance . , ■ • 


The testing of vehicle models which have been scaled down to 
some convenient size and operated in a simulated flight environment 
is one 'of the technicues which has been employed in ‘the acquisition of 
performance data for studying the base heating problem. 

Certain limitations e>:ist in scale modeling such as the size and 
operating capabilities of available facilities and the degree to which the 
fvill-scale vehicle can be modeled. The limnltations are compounded by 
the lack of unde r standing of the base heating problem which results 
in the inability to fully define all the necessary scaling laws. Two 
different test techniques were employed for the Saturn scale model test 
program. These were designated as the long-duration and short-duration 
techniques. The long-duration, technique consisted of mounting a scaled 
model of the vehicle in a wind tunnel test section and operating the 
nnodel at some simulated flight condition. In order to carry out such 
m-odel experiments on a continuous (or semi -continuous) flow basis,- 
complicated and difficult engineering problems encountered in the 
fueling, combustion, cooling, ignition, etc. of small rocket engines 
had to be solved. Also, maintaining a continuous high altitude 


Lj 


purripir.;; capaoiluy. ^ nis tecrimcue was lounc to oe too com 
costly to be consicerec for normal laboratory testing of Satt: 

Because of the difficulties with a continuous flow rocket test, a short- 
. duration technique was developed by Cornell Aeronautical Laboratory 
(CAL). The short-duration technique consists of mounting only a scale 
model of the base region in a vacuum chamber, and testing at'a simu- 
lated altitude. The actv-al test time for this technique is in the order 
of i few milliseconds, compared to the several' seconds of the long- 
curaticn tests, based upon results which show, that the base heating 
rates and pressures reach a quasi steady-state condition in less than 
a millisecond. Steady combustion is maintained during the time 
required for the expansion waves created by the rupture of a diaphragm 
to travel the length of the supply tubes and return to the' combustion 
chamber. 

One of the advantages of the short-duration technique is the 
nearly instantaneous establishment of a hot, steady flow, lasting 
for a few milliseconds, which has permitted the use of simple models 
without' cooling. 

Overall, the short-duration techniques appears to give the same 
results but at a lower price. Saturn I long-duration tests cost from 
'$50 to $150 thousand each (Ref. 5-1), whereas short-duration 
tests cost approximately 10 ne 'cent of this amount. 
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:ALING COhiSIDERATIONS 

To illustrate model scaling, consider some arbitrary differ- 
ential area in the base region s was showm in Section 2. 1, Figure 2-1. 
An energy balance shows that the total heat flux (qr) must equal the 
convective flux (q^) plus the absorbed radiant flux i.e.. 
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(5-1) 


°-t ” ‘ 

The cor;.vectiv£ and radiant flux can likewise be written as 


q_ = h(Tr - Tw) 


(5-2) 


and 


"^hlr ~ J ^p^-dAj — dAp 

A^ 


SpCrT^' dP^Aj —dAp (5-3) 


Ap 


The definition of these terms can be found in Sections 2. 2 and 2. 5, 
Equations 2-5 and 2-16, respectively. 

For the sake of brevity Equation 5-3 has been simplified and 
• written as if the plume radiation were a surface phenomenon. In 
some cases, the radiation can be treated as a surface phenomenon; 
however, in general, it is a. gaseous problem and must be treated 


as such. 


In model heat transfer testing, the obvious objective is to 
reduce the ph/sical diimensions of the prototype to some size which 
call be conveniently- tested. A relationship between the model and 
prototype can be established to relate the model heat transfer, to that 
of the prototype. However, with heat transfer problems as complex 
as base heating, the scaling becoTnes exceedingly complicated. To 
identify some of the problems which exist, the following is a discussion 
of factors which must be considered. 
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5.2. 1 Convective Modeline 


in scale rncceling the base heat transfer problenn, it is usual , 
to make the idealizing assumption that the model is an identical 
scaled version of the prototype and differs only in physical size. . 
Furthe rnnore, the convective heat transfer in the base region of both 
the model and prototype must satisfy- the well-known Nusselt equation 
(discussed in -Section 2. 2). 


m . n ' 

^nu = C Nr-e 


're 


(5-4) 


where 


N 


nu 


~ Nusselt Number - hx/k 


Nre 


- .Reynolds Number - pux/p 


x^randti Number = p.Cp/k 
C, rn, n - Empirical constants. 


The .ideal model assumption leads to the further assumption 
that the base gas density (p), velocity (u) and transport properties 
(k, p, are also identical or very similar. Thus wmen the model- 

and flight vehicle convective heating rates are compared using 
Equation a relationship is established, giving 


^m.odel / flight 

\ model 


. 1 -m 


h . 

-uiighr 


(5-5) 


which shows that the ratio of tiie heat transfer coefficients from model 
and flight vehicle are inversely- related to the ratio of their respective 
characteristic lengths (x) to the i-m power. Now, from Newton's 
law of cooiina, i. e. , 




u 


n I 


n ' 


U : 


u ■ 


n - 


i ^ 




i » 


( — ^ 


I ! 


! — ) f 


U : 


n 


V i 


rm 


LJ' il 


U ; 




t ! 


i 


\ 


uJ ’ ! 

I 


from Newton's 




an additional relation follows. 


oince 


h is not directly measured 


it is convenient to express r^quation 5-5 in terms which are measurable. 
Equation 5-5, using Newton’s law of cooling, becomes 



^-model 
. ^flight 




light 


-m 


-•^model 


(5-6) 


It has been assumed that the temperature differences for model and 
flight are equal, .i.e., ' 


<Tr - = (T, 


“ w 


^flight 


When dealing with convective modeling in. the past, it was 
assumed that the base region could be approximated by a flat plate 
(Refs. 5-2 and 5-3)* Under this assumption, the Colburn equation 
(Equation Z-6) for turbulent flow over a flat plate is used (Ref. 5-4). 
This equation has the Reynolds number raised to the 4/5 power ' 
(i.e., m - 4/5). Equation 5-6 then becomes 


/ -^model \ 


(5.-7). 


where x t/x.-._- is the scale of the model used in the test, 
model xiignt 

As an example, a summary of the models used in the Saturn 
program is shown in Table 5-i. By using Equation 5-7 and the scale 
given for each model,- the foliov/ing scaling factors result: 
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rcsui*:, when the prototype is scaled, rhe optical path length is also 
scaled, leading to a .change in tne ennissive power and therefore the 
incident radiant flux,' c-,,. ‘The latter assuines that the plume temperature 
profile is not altered in the- scaling process. As opposed to this the 
exhaust olume of a rocket engine operating at a low altitude has an 
afterburning mantle which affects the radiation. Because liquid 
rocket engines operate^most efficiently '-when the propellant ratio is 
fuel- rich, the exhaust products contain unbu’rned fuel. At some 
distz.nce downstream of the engine exit, the exhaust products have 
■ sufficiently mixed with the aiTibient air such that the unburned fuel 
ignites and forms the afterburning mantle (see Figure 2-1). ’The 
- distance downstream required to mi:-: the exhaust products and ambient 
air appears to depend upon the size of the engine, but this is not a linear 
rela':ion. As a result the afterburning characteristics of a prototype 
engine and a scaled engine appear to be quite different. This again 
leads to differences in the plume emissive power for the scale model 
and the prototype. In addition to the above, other factors such as carbon 
particle generation associated with tifd^^'iydrocarbon propellants, the 
engine and plumLC chemical kinetics,- and plume boundary layer develop- 
ment enter into the radiation problem and could be expected to affect 
'che scaling. In conclusion the relationship between the radiation 
associated with a scale model and that of the prototype would appear to y 
be difficult to establish at present. 

5.2.3 Gas Temperature Modeling 

• 3ecaus;e*of instrumentation difficulties, the recovery temperature 
on’ the rniodei cannot be determined directly. Therefore, some techniques 
have been developed that give an indirect method of determining the 
gas recovery temperature. One example of these techniques is the 
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‘'iriot-iiase iecnnicue (see Xex. o-d). . i-iiis particular recimique is . 
based unon Newton's law of cooling, which is • • . ‘ , 


^ conv. 


d( 


i- w ) • 


Trie basic premise c£ the' technique is that h is constant, that 
is, independent of temperature. Under this assumption a direct means 


of deter inining becomes evident^ i. e. , by varying the 


base teznpe ratur e, over a wide range in order to produce a significant 


igc in the measured q^onv * Now, since the gas recovery 


temperature is identical to the adiabatic wail temperature, can . 

be found by extrapolating the test data to a value of aero 


Actually, his a function of the temperature; therefore, a lower 


value of T-. would be indicated (Tv ). This technique,' however, 
^ "act 


dees give a good indication of the gas. temperature and can be used to 


. nci cate trends and possibly mag nitud e s . 


5.3 SOME APPARENT IMITATIONS IN SCALE MODELING -■ 

In practice, realistic limitations exist insofar as what can be 
scaled and what is beneficial to scale. Some examples of these 
limi ta tion s wi 1 h no w b e di s cu s s e d . 

Certain functional parts and components of the engine and 
vehicle are not duplicated. Examples would be the wall thickness of , 
the rocket engines which, if scaled accurately, would not contain the 
flow passages inside the 'wall. Therefore, the engines were generally 
larger than the ideal scale. 



Base gas temperature is related directly to nozzle boundary 
layer v/hich, in turn, is related to the engine size, contour,, and 
operating characteristics. 

Turbine cxlnaust gas discharge varies between model and. 
prototype. For example, ethylene turbine exhaust would more nearly 
simulate the actual gas condition on the prototype; hovvever, to obtain 
r:. steady level of turbine exhaust flow during combustion, the turbine 
exhau-st would have to be initiated prior to ignition. Also, for a 
few low-altitude runs where ethylene was used it was found to burn 
and coat the heating gages with soot. For these reasons, hydrogen 
gas was preferred over ethylene as the turbine exhaust gas. 

Regenerative cooling of 'nozzle wall is not simulated in the 
scale model engines . ■ This is expected to be a serious limitation 
in the case of the upper, stage engines (i.e., J-2 engines), where the 
effects of condensation on the inside nozzle wall can greatly influence 
the amount of energy which is reversed into the base region. 

Additional discrepancies which exist between scale model/ 
and prototype include the use of gaseous propellants as opposed to 
liquids, combustion chaimber injector design,' and the omission of 
an external flov/ in some of the booster stage runs. 

It appears that with the present state of. the art of scale rnodel 
testing, one model which successively reproduces, or models, 
both the convective and radiation heating phenomenas is not available. 
One ca'A be simulated only at the expense of the other. Therefore, 
in modeling vehicles, two different criteria need to be considered: 
one for radiation and the other for convection. 










The ability to correctly scale model test data to the full-scale 
vehicle requires au understanding of the heat transfer mechanisms 
in the base region as mentioned in the previous section. It has been 
assumed that the heat transfer ban be correlated with an equation of 
the Coiburn-Nusselt type; 


^ n m 
h.’nu = C jN:>^ hip 2; 


For similar surface flow proces ses between rocket nozzles and ba.se 
region, it can be shown that a functional relationship of the forms 

°base “ cham 

and (5-8) 

°^base ^ ^base 

exists and naight be used to predict the base convective heat flux. 

The major disadvantage of this naethod of correlating data from 
scale naodei tests is that it is only true for 'she same engine shape, and 
configuration. An example where the proportionality has been used in 
scaling thermal design criteria was in the base regions of the S-I, 'and 
S-IV vehicles. The K-1 engines used- on the vehicles went through a 
deveiopiaaent stage that increased its thrust from 165-, 000 Ibf on the 
S-I, Block I, to a thrust of 200, 000 ibf 01 a the S-IB. 

From data taken by the Cornell Aeronautical Laboratory , (CAL) 
(Ref. 5-6) on a four- engine configuration, the correlation obtained 
using Equation 5-8 was as shown in Figures 5-1 and 5-2. From 


5-i,' cxrronont (n) of Zquatior. 5-S is seen to vary bct\v eon 
0,75 to 0 . 9 s. T':'-c con'in'on value of the exponent associated with turbu- 

lent hear transfer over, a flat plate is 0. 8. Figure 5-2 shows the second 
correlation in Equation 5-S, The apparent power of 1. 02 is greater 
than the 0. 8 power which is -well established for flat plate turbulent 
flow. ■ . ■ 
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Therefore,, it seen that Equation 5-6 does not j.ccurately 
represent rhe base region of the vehicle. This is believed to be . 
caused by the assumption tha.t the recovery temperature is independent 
of chamber pressure in the development of these relationships. The • 
recovery temperature as a function of chamber pressure ior the same 
data, represented in Figures 5-1 and 5-Z is shown in Figure 5-3, It 
shows that there exists a relationship of the form T-p^ a where for 

“his particular case the average exponent is seen to be 0.20. Thus,* 
the assumption that recovery temperature was constant in developing ■ 
Equation 5-S was apparently incorrect. 

In light of these discrepancies, the model data was replotted 
using the Nus selt- Colburn type equation with Pr = 1 . These results 
are sho'wn in Figure 5-4. Figure 5-4 indicates a turbulent variation 
) in heat traiisfer with Reynolds number over a major portion 
of the base region. However, the, scatter in the data at the lower values 
of indicates a possible correlation v;ith laninar flow theory where 

^nu O' A sinailar anomaly ’was observed on the five-engine 

S-il model test. . ' • 

As an attempt to better understand the flow phenomena occur- 
ring during full flow reversal, some very detailed, model studies . 
have, been conducted. The results of one such investigation are 
reported in Reference 5-7. During these tests probes were inserted 
into the base region- of a four-engine coruiguration. These probes 
measured "he static and impact pressure as well as velocities. The 
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velocities were measured using hot wire anemometers. Some typica.i 
results d£ tViese tests are shewn in Figures 5-5' through 5-7, These 
figures shov.' the static and impact pressures and the resulting Mach 
numiber obtained bv combining these results. 




Several interesting observations can be made in studying these 
figures* Figure ^5-7 shows that supersonic flow can apparently exist in 
the base region. The static pressure shown in Figure 5-5 shows a sharp 
increase as the flow approaches the heat shield (e.g. , compa.re the 0. 5- 
inch curve with the 1. 0-inch curve in Figure ’5-5). This trend, if not 
the iTiagnitude, is in agreement with what would be anticipated, with 
supersonic flow conditions {i. e. , a standing shock above the heat shield). 
Measurements taken along the base plate also indicated sonic velocities 
as the flov/ approached the vent area. Efforts to show the existence of 
the base plate , standing ‘Shock using a laser-Doppler technique (Ref. 

5-S), however, did not provide concrete results to confirm the existence 


of the shock. ... 

Another, interesting .factor ‘which can be obtained from these 
results is the following. Isentropic flow is normally assumed in the 
base region, and isentropic flow relations are used in the analytical 
models. In isentropic flow, no change in total pressure occurs unless 
the flow is supersonic and moves through a shock wave. As the 
flow :r-oves through the shock, the entropy changes.- In checking 
the data shown in Figures 5-5 through 5-7, one would expect to find 
the entropy changes indicated by the change in impact pressure to 
correspond to the entropy change across a normal shock. What is 
found, however, are entropy changes which are 3 to 10 times, greater 
than the entropy change across a normal shock. This could mean 
that the'flow in the base region has much larger viscous losses 
associated with it than has previously been assumed. 
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encrine S-II scale model test 


However, from studies on four- and five- engine configurations (Reis. 


5-7 and .5-9), it is interesting to note the dissimilarities between the 
two* . Figure 5-8 shows a graphical representation of the velocity 
profiles which appear to exist. The t^mical velocity profile of the 
four-engine cluster appears to resemble that' of a jet impinging on a 
flat plate; however, the five-engine case shows that a inajor portion of 


the flow appears to leave the base region before stagnation on the base 
niatc -occurs, ■ 


5,5 MODEL- DATA AND FLIGHT DATA COMPARISON 

5.5.1 S-IC Stage Model Test Results 

A model test program was initiated early in 1962 to aid in ■ 
predicting the base heating environment on the Saturn V/S-IC stage, 

A short-duration (shock tube) test technique was chosen because it 
provdded a much sizTipler and cheaper method of testing. Confidence . 
in the technique had been obtained earlier by the favorable comparison 
between the long and short duration test results. A l/45th scale' model 
of the S-IC stage base region was designed by Cornell Aeronautical 
Laboratory (CAL). 

Tests were conducted in CAL’s high-altitude chamber at altitudes 
between 125, 000 and 205, 000 feet with no external flow and in the CAL 
transonic tunnel between Mach 0.6 and 1.2. The NASA/Lewis Research 
Center transonic and supersonic tunnels were used to obtain model 
data in the Mach 0.4 to 3.5 range. 
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• The model used GOX/ethyiene to simulate LOX/RP-1 pi*o- 
peilants. Hydrogen was injected into the main exhaust flow inside 
each of the nozzles at the 10: i area ratio point to simulate, the turbine 
exhaust. 

Thin-film heat transfer gages (calorimeters) were employed 
for the measurement of short-duration heating. Because of the short- 
time duration of the model test, the gages were recording a cold w?Hl 
heating rate. 

As shown in Figure 5-9, the scale model radiation heating 
rates were negligible as compared to those measured during the 
flight tests. This difference .is probably due to the lack of particulate 
carbon (soot), the smaller radiating gas volume, and the afterburning 
mantle location that is associated with the scale model engine plumes. 

Model data total heating is. showii in Figure 5 .-10 for the heat 
shield. It is made up almost totally of convective heating and follows 
the trend of the flight data. The model convective heating presented 
in Figure. 5~11 is observed to be much greater than the flight data. 

The scaling factor for convective heating discussed in the previous 
section was applied to the model data and is also shown in Figure 5-11. 
The model data did not indi.cate the convective cooling which took 
place on the flight vehicle early in flight; however, the tests did indicaie 
the trends (increases and decreases) in the base region forTater times 
in flight. For comparison purposes the convective heating values 
obtained from the aspiration model and the flow reversal nrodei 
developed at Teiedyne Brown Engineering are presented. 

Base pressure measurements from model and flight are com- 
pared in Figure 5-12. The shapes of the curves are very similar; 
however, the model base pressures have not been corrected for the 
differences in efficiency between model and flight engines, scoops, 
and flow deflec'cors. . • ' 
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ZTiigrit and model gas temperatures are compared in Figure 5-13. 
The flight data' shows large gradients in the gas temperature near the 
heat shield. The variations in the model data are significant, however, 
it is evident that the flight and m.odel data have comparable magnitudes 
of temperature at the higher altitudes. Data for the lower altitudes 
(or :vlach numbers) were not reported, since CAD felt that the data 


obtained for tin 


:onditions were not reliable. 


S-IB Stage Model Test Results 

A model test program on a 5. 47 percent scale of the S-IB stage 
booster v/as conducted at Arnold Engineering Development Center 
(AEDC) to determine the effect of turbine exhaust gases on base heating 
Short-duration- techniques were used in the test program.- Tests were 
conducted at >4ach numbers of 0, 8, 1. IS, and 1.63; however, because 
of tunnel flow disturbances caused by the mass flow addition of the 


roc-kets into- the 


tunnel, the data at Mach 0. S was believed unreliable 


and omitted. Therefore, only data for Mach 1. IS (30, 000 feet) and 
1.63 (40, 000 feet) is available. COX/ethylene was used to, produce the 
main engine exhaust, and hydrogen gas at 100® F was used to simulate 
the fuel-rich turbine exhaust gases. 

On the. S-IB stage there are two areas 'of importance in studying 


cne case 


vi r o nm e nt : th e h 


hield and the flame shield. 


For. this reason the comparison will be broken into these two categories. 

The same difficulties mentioned in the preceding section on 
the S-IC stage, in simulating the base heating, also apply 'here. 

The total heating rate and radiation heating on the heat shield *. 
are presented in Figures 5-14 and 5-15» The heat shield thermal 
environrneht of the model was relatively unaffected by the turbine 
exhaust configuration (i.e., whether the turbine exhaust is exhausted 













ii four 'inboard eni?ines, or exhausted outboard through the stub 


j , 


Figure 5-16 shews the convective heating obtained fronr 
■'igures 5-14 and 5-15 for the heat shield. The scaling factor was 
to the model data and this is also oresented in Figure 5-16. 


ex ',J u 11. kS \ 


in r Igures o- 


T7 and 5-18 the comparison ,'beuveen the flight, 
and naodci total and radiation heating data is'presented. As illustrated 
in Figure 5-17, the simulated turbine exhaust gas has a definite effect 
on the lotal heating value. The hydrogen gas' used to simulate the 
turbine exhaust was at 100° F, whereas the flight vehicle’s turbine 
exhaust is approximately 1, 500° F. Therefore, in the case when the 
turbine exhaust is exhausted around the. flame shield, the gas has a 
cocling effect in the model where it may actually have a heating effect 


in flight. 


. No appreciable difference was found in the model’s radiation 
heating, with or without simulated turbine exhaust gases, or because of 
different configurations used to dump the turbine exhaust gases' (Ref. 
5 - 3 ). * ■ . 


5.5.3 S-II Stage Model Test Results 


Parametric model tests were conducted by CAP in 1962 using a 
1/25 scale model of the S-II stage and short-duration techniques. The 


niodel made, provisions for variation in engine pitch circle, heat shield 
location, and interstage skirt location. 


The model v/as moved from CAL to the larger altitude chamber 
ar ‘viarshaii Space Flight Center (MSFC).- This move was required in 
order that the thrust structure heating rates could be better analyzed. 
In the CAL test, shock reflections from the tank walls disturbed the 


flow within 3.4 milliseconds. 
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n this tin^e interval the rhrus' 
~eec*oa*ce ci.irie uo SLc:.oij.i-i-c* 
:in“ios in the order o: 6 to 7 m 


structure heating rates did not have 
h'ISFC, the larger tank provided test 
lliseconds before flow disturbance. 


The test on the S-II stage was limited to the capabilities of the 
chambers. The CAL chamber has a vacuum capacity of approximately 
b nmlcrons of mercury (Zop, 500 feet) which is only 170 seconds into 


: C 'f' C: ■'r 


jents 10 to 20 seconds of S-II stage flight time. 


Lccause of the limiitations imposed by these facilities, a simulated 
altitude of only 266, 500 feet was obtainable. 


The comparison of base heating data for the S-II stage can be 
presented as follows: 1) Heat shield thernral environment, and 

2) Thrust cone thermal environi-nent. Figures 5-l9 ^nd 5-20 compare 
flight and model data on the heat shield, and Figures- 5-21 and 5-22 
cemnoare data on the thrust cone. As shown in these figures the heat 
shield experiences a more severe heating environment than the thrust 
cone. Looking, at the flight data, the thrust cone is seen to have a 
negligible convective heating rate during flight (approximately 0. 02 to 
0. 04 Btu/ft“-sec). 


i^fure o 




3 presents a comparison of the flight and model 


convective heating rates for the heat shield. The model data has been 
scaled assuming a fiat plate correlation. As shown, once the model 
data is scaled, it is in better agreenrent with the flight data. The 

also preserxted for . comparison. 


analytical value ! 

Some typical test data for the S-Ii scale model obtained by 
CA.L are presented in. Figure 5-24. As noted,- the values for heat shield 
heat transfer and pressure level off after about 0.5 millisecond. The' 
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GURE 5-19. S-II STAGE FLIGHT AMD MODEL HEAT SHIELD TOTAL HEAT FLUX 
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■ FIGURE 5-24. TYPICAL S-II STAGE TEST DATA • 
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thrust cone heating rates on the other hand do not appear to re?,ch 
a steady- siate condition during the test. This is thought to be caused 
by the size limitation of the chamber. The test apparatus was moved 
to NASA/hlSFC to investigate the thrust. cone in more detail. Results 
of this test compared to CALL'S data are shown in Figure 5-25. The 
CAR data gives the impression that it might rise to a higher heating 
level if it were not interrupted by tine returning shock. The fact that 
heatimr rates rise at a slower rate in the- larger tank makes it apoear 
that the heating rate may be caused by a flow phenomenon characteristi 
of the chamber rather than the effect of exhaust gas reversal from the 
model. 


5*6 CONCAUSIORS REGARDING MODEA TEST RESULTS 

Some conclusions of the scale Aiodei testing are listed below. 

o In general, measurements obtained directly from 
tests performed on the model appear to indicate 
a higher base thermal environment than that 
measui-ed during the flight. 

o o 

o The .application of scaling factors to the convective ■ 
flux gives better agreement between the nnodel and 
the full-scale flight results. Thus> it appears 
that the model results should be scaled; however, the 
exact scaling laws that apply are unknown and are 
probably more complex. than those which have been 
used. 

o Judicious interpretation of the results- of model testing 
appears to be required, as evidenced by the wide’ 
scattering of test data and by the methods used in 
certain tests (e*g*, simulation of the turbine 
exhaust gases of the S-IB flame shield). 

o In most cases, the trends with respect to con *7 • 

figuration changes appear to be indicated correctly. 
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c V/ith respect to the Saturn model test, 
points were nored: 


following 


A The base pressures were predicted fairly 

accurately and correlations have been developed 
using model and flight data in other references’ 

(Ref. 5“ 11). 

A The convective, cooling that was evident early 
in the.-S-IC flights- was not indicated by the 
model. 

A Model heating rates that were obtained for the 
upper, stages were in better agreement with the 
flight -re suits; however, the gas temperatures 
'that were indicated by the hot base technique 
were significantly above those indicated during 
the flight tests . .. . • 

c Seme of the problems and inaccuracies in scale model 
testing appear to be caused by limitations in'the 
facilities and the lack of s.imulacioh of certain 
functional details of the vehicle, rather than by the 
test technique itself. 

b A problem appears to exist in translating scale mxodei 
test data into effective flight vehicle design criteria. 

In ‘order topredict vehicle surface temperatures, it 
is necessary to know the convective heat transfer 
coefficient, the recovery temperature, and the inci- 
dent radiation. The exception to this is when the 
total heat flux is of such a magnitude that the surface 
reradiation will cause equilibrium to exist at a 
terxiperature below the limit of the material. The 
latter was the situation for most of the Saturn upper 
stages. For all of the booster stages and for many 
specific problems on upper stages, it was necessary 
to estiimate the probable gas temperature and- heat 
transfer coefficients. The desired design criteria 
could be established much more conclusively if these 
parameters could be obtained .from the scale model 
test results. Therefore, the- following is recommended 



LJ 



K 


Mere cf-ort should be devoted to developing 
techniques for measuring base recovery 
rc'mperatur es and evaluating convective heat ■ 
transfer -coefficients from scale model test 


ivlore basic studies such as thos e repres ented • 
by Fceferences 5-7 and 5-9 should be carried 
out to investigate basic parameters such as 
base flow fields, velocities, recovery tem}:)era.- 
turc,. and convective heat transfer coefficients. 



A more thorough exp er inrental investigation of 
the scale model to flight vehicle scaling lavvs 
should be pursued. 
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. 6. BASE KEAT TRANSFER DIAGNOSTIC' AND 
PREDICTION TECHNIQUES 

6. 1 ASPIRATION PHASE MODEL 

The Saturn flight test results have indicated that aspiration of 
atnbient air into the base re9;icn exists 'durins; the first part of the 
flight* During this phase, airabient air flows into the base region. If 
the radiant flux is high, the relatively low temperature ambient air 
provides convective cooling. If should be noted' that the Saturn I and 
S-IB flame shield region was an exception to the above. A semnl- 
empirical analytical model which approximates the base convective • 
environment during this phase has been devised and is described below. 


In Section 2 it was shown that in order to predict the convective 
flu:-: the following para-meter must be determined: 


q^(r) = f[li(T), T^(t)] (6-1) 

h(r) = f[c, m, ,n, Pq, ^1' , ’ ' (6-2) 

T^(t.) = £[To(r)] .,(6-3) 

where - 

- convective .flux 

h - convective heat transfer, coefficient' 

T.^ - gas recovery temperature 

pQ - stagnation pr*essure 

?.c - static -pressure 


tac-nation tem*oeraiure 


c^*.aracteristic dimension 


' > * * * ; - * 


- e:-:nerimentailv determined coefficients. 


Figure 6~1 is a' simplified schematic of the proposed model. In 
this .modeij the flow is shown to make a 90-degree turn at the edge of the 
base region and flow radially into the center. In reality, the flow would 
become detached at the edge and probably vortices would be shed into the 
base region. This flow model v/as ■ selected partially for the following 
reasons: 

o The sui^ersonic engine exhaust jets produce a . 
definite aspiration effect upon the rocket base, 
tending to draw the ambient air into the region. 

o The use of peripheral scoops pn some of the 

vehicles tended to divert air into the base region. 

o Prebabiy as a result of the above, measurements 
of rocket base heat transfer indicate an. inverse 
function of the base radius. , This can best be * 
explained by a flow, which is radially inward 
across the base. 

Using this simplified model the parameters of Equations 6~1 
through 6-3 can be approximated as follows. The trajectory of the 
vehicle is always known; therefore, the freestream ambient conditions 
ar-e also imowrx. The 'most utilitarian model of the base aspiration 
phase heat transfer would be one which relates the base heating to the 
freestream parameters. ’ ' ■ ■ . 

Figure 6-2 shows the experimentally determined pressure in the 
separated flow region of a bluff body and a. .back step taken from Ref- 
erences 6-1 and 6-2. Figure 6-3 shows ho v/ the bluff body base pres- 
sure agreed with the corrected -pressure measured in the base region 









REESTREAM MAGH NUMBER 

IMENTAL BASE TO 'AMBIENT PRESSURE RATIO 
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Cl z.ic j^arurn 


AS-dOI flight. 3y letting the ambient and bluff body ba: 


pressure r.cpresent, respectively/ the stagnation and static pressure of 
the radial component of flow into the base region, a base flow Mach 
number can be commuted from 


M/(r) = 


) . 


(6-4) 


n • 


Li ’ f 


2v!-j^(t) - Mach number of the radial inward base flow at 


i M 


- center base pressure at time r. 


n i 


PUt) ' - local freestream ambient pressure, at tim*e t 


To assist in evaluating the base-to-am.bient pressure ratio, the bluff 
body base pressure curve shown in Figure 6-2 was curve-fit between 
freestream Mach numibers 0 = M^^ =1.5 to obtain 


1 - 0. 13.5 


0 S Moo^ 1.5 


(6-5) 


'igure 6-4 shows the corn-parison between the base flow Mach 


i I i 


nurnber obtained using Equations 6-4 and 6-5. and that evaluated from 
the Saturn AS-501 and A.S-502 measured -pressures. ' The figure illus- 


1 I 

lJ ij 


rated that the computed Mach number agrees reasonably well with the 
fPtained from the measured data during the Aspiration phase. 


Based on. Saturn S-IC flight results, the stagnation temperature 
me base region is assumed to be the same as the freestream; thus, . 


n ' j 
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= '^cc. 
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T:ic base recovery temperature can be obtained from 


-00 i - 


, %/ - i . . 


-'Vj 


(&-7) 


Sogin (Pv.ef. 6-3) found that a recovery factor, r^ = 0. i, gave good 
agz'cc'ment with heat transfer measured in base regions of bluff bodies* 
Thus^ using = 0* 1 and v = 1.4, the recovery temperature can be 
obtained from 


^ (1 -t O.OZ Mod ) 




(6-S) 


"Wlicn evaluating transport properties, the static temperature of the 
base gas c 2 in be obtained from 


/ 


-s = 


/, -.xm vrr 


(6-9) 


Combinina- E-cuations 6-9 and 6-6 gives 
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I v=_l 


M-k“ 


(6-10) 




oort properties for air should be evaluated at. the film temper a- 


ire given by 


■-i-'S “ -w- 


(6-11) 
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suriace tem'oerarare . 


r:ased upon the simplified flow model of Figure 6-1, the characteristic 
dimension is taken as 


X = - r . 


( 6 - 12 ) 


Using the. equation derived in Section 2, i. e. , 


= h(T,- - T^) . 


(2-5) 


v/nere 


/ 


h = C 


n 1 — m m 

Cp/ Kf- --p‘“u"‘ 




■ !^f 


m-n 


j 


X 


(2-7) 


and coefficients' m = 0. 8; n = 0. 333; the convective flux to the S-IC heat 
shield calorimeters C25-106 and C26-106 was computed. The computed 
results, compared to the mieasured convective flux, are shown in 
Figures 6-5 and 6-6. The coefficient C = 0.03 is equivalent to flat 
plate, parallel turbulent flow: heat transfer. For the aspiration phase 
convective cooling, using C = 0. 03 would give conservative design 
results. Based on the results shown in. Figures 4-28 and 4-29 and 
Table 4—1 for the heat transfer in separated regions, -a higher coeffi- 
cient may be justified. Therefore, the convective heat transfer was com.- 
puted using coefficients C = 0. 06 and C =,0.09. These results are also 
shown in Figures 6 -.5 and 6-6.. 
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The figures illustrate that the convective Ulo: computed using the 
-tion phase analytical model is in reasonably good agreement with 
ensured flux evaluated for calorimeters C25-106 and C26-106; 


I \ 1 . 


.C25-106 thjan for C26-106. The figures also shew that the asoiration 
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nodel will yield 
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PLUMi-J AND AFTERBURNING EFFECTS 


In studying booster base heating and related problems, the 
.heroes of the engine exhaust olume s. beconre invoortant for several 


xts.ons. -Such factors 


as the rime 


ind/or altitude at which' flow rever- 


• sai. srarts the apparent coincidence of increases in radiazit flux 
during this period are related to the plume shapes and interactions . 

To further understand these relations, a cursory study of actuahplume 
shapes, predicted inviscid shapes, and adjacent plume interactions, was 
initiatc.d with the following preliminary results. 

At present':, sevora.1 methods of predicting inviscid plume shapes 
are available. Throughout this study an -approximate method developed 
by Charwat (Ref. 6-4) has been used. Comparisons have indicated that 
this approxinnaticn gives good results compared to the more involved 
method of characteristics . 


rr ; 


Figure 6-7 sriows a dimensionless plot of a single plume bound- 


ary wnich Wc 


Lculated using the method nnentioned above for an F-1 


engine with a 3:1 area ratio. Two different types of input were used, in 
the calculations. The first utilized the NA-SA/Lewis Thermochemicai 
Fro ‘u ram. to evaluate the nozzle exit condition's, whereas the second used 
the one -cimensioriai isentrcpic flow theory to obtain the exit conditions, . 
From photographs taken during static tests of this same engine, the 


o 




actual plume boundary has been scaled.and nondimensionalisod and is 
shown in Pi^'ure 6-7 alon<^ with the inviscid olume shapes. It can bo 
seen in the figure that the inviscid plume boundary predictions are sub- 
stantially less than the real plume boundary, especially in the after- 
burnin{r mixinc^ zone. 


To show how this affects the base heat transfer, 


the gas temper a- 


iures and pressures measured in the flame shield and heat shield 
regions during one of the Saturn I flight tests are shown in Figures 6.-S 
and 6-9. As shown in Figure 6-8, flame shield flow reversal appears • 
to start at approximately 25 seconds, whereas the heat shield flow reversals 
start at approximately 60 seconds. (Note. ■ These resvilts have not 
been normialized to AS -501 trajectory.) Figure 6-9 shows the corre- 
sponding pressures to be approximately 10.2 and 3.4 psia,. respectively. ■ 
Figure 6-10 shows the dimensionless plume shapes of the Saturn engines 


predicted using Charwat's method fer various pressures. The center- 
lines betv/een adjacent engines are also shown in the -figure. 'By studying 


Figure 6-10, it can be seen that the inviscid plume prediction would not' 
have predicted plume interaction corresponding to the flame and heat 
shield pressures of 10.2 and 3,4 psia, respectively. Apparently, the 
tnick nm'xing and afterburning layers are interacting and reversing 
much earlier in flight than would be indicated by the inviscid plume 
predictions. Therefore, using plume shapes which have been generated 
using inviscid flow theories to predict plume impingement may not be • 


r epre s entative of the real c:ondition. Further analytical work to 
Incorporate the mixing zone in the prediction of plume shapes is 
indiciated here. 

A simkiar analysis can be conducted on the F- I (Saturn V, S-IC) 


engine. The anaiyuical plume shapes for several altitudes, for the F-1, 


are presented in Figure 6-11. 
























Single engine plume shapes for the J-2 (S-II stage) engine have 
been calculated and are presented in' Figure 6-12. Afterburning is not 
a problem on the .upper stages, as it. is- for boosters. Therefore, it 
should be mentioned that the inviscid plume shapes should be more valid 
for upper stages. 


6.3 BASS RADIATION 

As -discussed in Section 2, rocket exhaust plumes are noniso- 
thermal, nonisobaric, and noniiomogeneous gas bodies. The equation 
w'nich governs the radiation from such bodies is 


Cir = \ \ P K,^ I-bv exp (- \ p Ky dS') 


X cos dS duJ dv . , 


(6-13) 


wner e 


rident radiant flux 


local densi'cy 


- spectral absorption coefficient 


- Planck 


lancK rune non 


- optical oatn 


cj - solid angl( 


V ■ - wave numoej 


- polar angle (Figure 2-: 


For conditions where the radiation can be considered a surface 
phenornex>on, the equation for radiation between a finite radiating plume 


o uriace e 


clei-ient and a differential receiver is 


6-19 









(6-14) 


^ir‘ ~ A Aq — AAr^ 


^v:ier e 


O' ■- . plume blackbody emAssive power 


• ^ AAS - AAp. 

6.3*i Gaseous Radiation 


- view factor between A Ac and AAv 


A number of methods have been' developed. for predicting gaseous 
radiatio'A and Reference 6-5 contains an excellent literature survey on 
this- subject. The most successful method involves the use of several 
computer programs and was developed for or by the MSFC Aero- 
Astrodynamics Laboratory over a period of years. Briefly, the 
calculations are carried out in the following way. 


The chenaical composition of the reacting propellants is obtained 
using a thermal-cheinical computer program similar to that described 
in Pv^eference 6-c. This program generates the chemical composition 
of the gases as an input to the method-of-characteristics flow program. 


A method-of-characteristics (Ref. 6-7) flow program in con- 
j“unction with the thermal- chemical program is used to generate the 
inviscid flow fields of the engine and plumes,. This program can also 
■generate the internal shock structure of the plumes to a limited degree. 
Since the method-of-characteristics is restricted to supersonic flows, 
the Riernann waves (normal shocks) which are known to exist at points 
along the plume axis cannot be generated with this program.. Also, no 
computer program is known to exist which can readily generate the 
three-dimensional shock waves which would exist as a result of the 
piume/plume inter actions . In some cases, approximations of the 
shock structure in tlaese regions have been used and patched into the 
inv'iscid supersonic flow fields with some success. 


} 


n 


Si:‘:Cc the me th od - c I - char actc r 1 s t i c s program output is the .. 
irer ptodyramic and flow properties at a point where right- and Icft- 
runntng characteristic lines cross, it is not usually suitable for direct 
input into the gaseous radiation programs. A, separate compute r .pro- 
gram' vailed JUGGLE (Ref. 6-S) is used to prepare the n*iethod-of- 
cr.aractcristics output for use in the radiation program. This program 
generates a systematic miap of plume thermodynamic propertie s .as a 
• function of axial distance downstream of the exit and radial distance 
■from the ttxis. This is accomplished by iterating between the four 
nearc.SL points of the* method -of- char acteristics output. 

If afterburning along the plume boundary is significant and is 
to be included in the radiation problem, the thermodynamic and chemical 
compositicri of the afterburning mantle must also be supplied. One of 
the many techniques w'hich have been deyeloped for this purpose is 
discussed in Reference '6-9. This technique uses the. streams tube 
approach, viscous mixing, and finite rate chemical reactions to com- 
pute the temperature and gas properties of the mixing boundary layer. 
This layer is then superimposed on the ihviscid flow field.- 

Otxcq the thermodynamic properties and chemical species are 
spatially defined, the thermal radiation can be computed.- The radiation 
calcuiariions are essentially the, numerical integration of Equation 6-13. 
Reference 6-iO describes a computer program which has been developed 
for this purpose. This program is a ,volu*imetric integration from an 
arbitrarily oriented differential area through'the plumes and surround- 
ing media along each line- of sight. The program uses band-averaged 
abs.orption coefficients (xy in Equation 6-13) and the Curtis-Godson • 
method of approximating nonhornegeneous gas properties. 

The techniques just outlined are theoretically sound and appear 
to be the most advanced if not the only complete technique presently 
available for computing plume radiation from basic knowledge. 


n 

LJ i 


I ^ 


L] ! 


u i 


r~ 


n ' 


U 


! 1 


i I 


t .1 ' 


J : 


( — ' 

U ^ 





Figure 6-13,. v/hich was extracted frona Reference 6-il, shows some 
-jynicai results which can be obiained -using this method to ’predict the 
“radiation from a single piu'cme. Given the complexities associated with 
the three- di'mensional flow field ai^d shock structure of multiplumes, 
chemical reactions in free bouerdaries and the unknown nature of particle 
•sices, concentrations, and densities^ it would appear that these tech- 
niques would not easily lend themselves to a simple solution to the prob- 
lem. Kov/ever, reasonably accurate predictions of the incident flux to a 
point in the base region can be obtained provided the . expenditure of the 
time and effort implied in using these tecliniQues is not prohibitive. 

As is sc often the case, a simple engineering miethod (it e. , one which 
is reasonably accurate most of the time and simple to use) would be 
very beneficial in solving this problem, 

6.3.2 Surfac e Radiation' Approximation ■ - 

For certain conditions , . it is possible to treat the plume as a • 
diffusely radiating surface. The physical implication of this approxi- 
nration is that the absorption coefficients associated with the gas 
.in the outer layers of the plume are so large that the gas is essentially 
opaque for any signiiicant thickness. This situation exists when the 
gas is heavily laden with particulate materials, ’ 

By treating the radiating plumes in this manner, the problem 
becomes mathematically simple as shown by Equation 6-14. The problem 
with/ this appr oxirnation is that the spatial variation in emissive power 
of the plumLC m.ust be known in order to predict the incident heat flux. 

Ofien rhis.can be deduced from radiation measurements taken during 
en^iine static test firings. 
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To aid in calculating the heat transfer between a diffusely 
plume and an arbitrarily located area in the base region, 
a computer program has been developed (Ref* 1-6). The details of 
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oi tp.is program arc given in the appendix. The program nas tne capa- 
bility to compute the radiation zr orn. multiple plunees, including any 
blockage by intervening plumes or nonparticipating surfaces-. 

To show how well the calculated radiation agrees with the 
measured radiation, comparisons were made with the Saturn S-IC 
flight test results and are shown in-Figures 6-14 through 6-21. The 
F-] engine plume emiissive power, shown in Figure 6-14, was obtained 
bv extr anolating measured results taken during static tests of an H-1 


(Saturn S-I) m 


power to tnc 


In extrapolating the H-1 engine plume emiissive 
ngine' plunae, it was assumed that the dimension- 


less distance, X/D, would be the samie for both engines. The only justi- 
fication for this is that both engines used the same propellants (LOX/ 
RP-1) and the afterburning appeared to start at approximiately X/D 1 


:or both. 


X It is anticipated that better agreeiment would have been 


obtained if the actual F-1 emiissive power-were available. The plume 
shape was obtained fromi photographs taken during the static tests. ■ 


ine plumes 


re divided into short cylindrical segments over which 


missive 


power could be assumed constant as shown by the dashed 


line in Figure 6-14. Since the available literature indicated two possible 
levels of the plume emissive power (Figure 6-14), both-of these were 


rwo heat shield locations corresponding to the heat flux meters, C60-106 
and C6 1-106 shown in Figure 6-15, and three engine locations corre- 
sponding to C57-101, C58-105, and Ci50-i01 shown in Figures 6-15 
and 6- 16. The calculated and mieasured results are shown, in Figures 
6-17 through 6-21. ■ For the .calculated results, the upper and lower parts 
of the band shown in the figures correspond to the upper and lower em.is- 


uve power cur 


ves of Figure 6-14. Because of la'onch pad' deflections of 


the F-1 engine plum^es, the ccimputed flux should be compared with the 
flight heat flux measured after lift-off (possibly 5 to 15 seconds). • 
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FIGURE 6-17. S-IC MEASURED AND CAL.CULATED SEA LEVEL RADIANT FLUX 
AT C60-106 









FIGURE 6-19. S-IC I'lEASLIRED A!!D CALCULATED SEA LEVEL RADIANT, FLUX 
TO Cl 50- 101 
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FIGURE 6-21. S-IC MEASURED AND CALCULATED SEA LEVEL RADIANT FLUX 
TO C57-101 









previously s:a::cd; the shortcomirig of tnis nietnod is its. 
reliance on che measured pluine emissive po\ver. Ideally the plume 
emissive ’power would be measured with a scanning radiometer from a 
position which offered complete visibility of the plume. However, ■ 
ihe accessibility of the plume and ihe environment often make this 
impossible. As a result, heat fi’ux meters located in close proximity 
to the engine exit and measuring only the irradiation to a point must be 
relied upon to estimate the emissive power. The emissive power is 
estiniatcd by iterating to determine what emissive power as a function 
distance will give agreement with the measured incident flux. 


6.3.3 Radiation at Altitude 

Because of the empirical nature of the surface radiation model 
and its dependence upon radiation measured during static testing of * 
a single engine ai sea level, it is not applicable to higher altitudes. 

One empirical technique which was used with reasonable success in 
predicting the Saturn booster base radiant flux is the so-called '‘radia- 
tion decay m,ethod’h This method consists of predicting the radiation 
at act altitude on the basis of previous flight test results. Some* of the • 
Imiitations associated with this technique have been discussed pre- 
viously in References 1-4 and 1-5. The general idea in using this 
technique is to select a vehicle with a configuration and propellants 
wkich are similar to those of the vehicle to be designed. Fortunately 
the inflight radiation characteristics in the base regions of enough 
dissimilar vehicles have been measured to make this possible ‘most of 
tne time. The base radiation characteristics, as a function of altitude 
of several vehicles, are shown in Figures 6-22 and 6-23. 

Because of the -similarity between some of the current space 
shuttle configurations and the Titan IIIC, it is anticipated that the base 
radiation characteristics will be similar to those shown for the Titan 
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me in. Figures 6-22 and d- 23. Although the Space Shuttle radiation is 
still being investigated, an indication of the level can be seen in the 
comparison between the' Saturn I, S-IC, and Titan IIIC' (Ref. -3-S) 
measured radiant flm: shown.-in Figure 6-24. 




6. 4- FMGIixE CRYOPRGPEFL-AiXT EFFECTS 

' ' With respect to the base heating problem, ' engine cryopr opcllant 
effects . are those effects which are caused by the extremely low- 
ternperature cryogenics. The 'most interesting example of this effect 
v/a's noted during, the vacuurn chamber test firings of the J-2 engine. 

Film taken during the tests indicated that a condensed liquid and ice 
* layer existed near the wall at the exit of the J-2 engine. At first it 
waa thought to be inconceivable that condensation could exit, con- 
sicioring the engine chamber temperature (and therefore the gas total 
temperature) of approxinrately 5, 700° F which is several thousand 
degrees above the temxperature at which condensation can be. formed. 
Kewever, considering the j-2 engine cooling circuit (Figure 6-25) in 
which liquid. hydrogen at a temperature of -423° F is forced through the 
wall of the engine, it was concluded that condensation was not only possi- 
ble but probable. In conjunction with this problem, a cursory analysis 
of rhe effect the condensation could have upon the S-II stage base heating 
was undertaken. The results of this analysis were documented in Ref- 
erence i-2. 

♦ Figure 6-26 shows how and where engine' condensation occurs 
and rhe internal wall temperature of the J-2 engine. At a point approxi- 
mately 40 inches downstream of the engine throat the wall temperature 
drops sharply. The Inydrog.en inlet to the cooling circuit and the turbine 
exhaust gas outlet are, both near this location and contribute to the sham 
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FIGURE- 6-25. THRUST CHAMBER COOLANT CIRCUIT 










is the satura- 


decrease. in v/ali temperature. Also shown in Figure 6-2 
'den temperature of the euhaust gases. Since the wall temperature is • 
below the saturation temperature of the gas, condensation can. occur. 

The study of the effect condensation could have on the S-II 
stage base heating was undertaken because, of the difference in gas 
temperature' predicted by the scale model test (~2500''F) and the in- ^ 
flight measured temperature ('^1000° F). Since the source of the 
reversed flow is a small portion of the fluid which'was originally in 
uhc nozzle boundary layer, it v/as postulated that if a condensed liquid 
layer existed it could revaporize upon leaving the nozzle and signifi- 
cantly reduce the base gas temperature. Other factors such as the 
energy loss to external engine surfaces and to the heat shield were 
also considered. The results of the study are indicated in the. follow- 
ing paragraph. 

The energy loss to me external surfaces of the engine could 
have reduced the m.ean temperature of the .base gas a maximum of 
apqjrcxirnately 200 ""F. The heat transfer to the heat shield could have 
reduced the nrean temperature of the base gas a maximum of 100® F. 

To reduce the mean temperature of the base gas an additional IZOO^'F 
would require a 2S to 43 percent moisture content in the reversed mass 
flew. Using a laminar flow theory (the only technique presently avail- 
able for computing condensation for the J-2 engine conditions), it was 
shov/n then more than the desired amount of -condens ation could have 
been produced in the J-2 engine, -As a result, it was concluded that 
engine wall condensation could have affected the S-II stage gas tempera- 
tures and heating rates. 

The possibility that the Space Shuttle engines might produce con- 
densation in amounts which could also affect the base heat transfer has 
been considered briefly. The latest Space Shuttle engine wail temperatu 
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j:rec.:.cuo:\s, obtained iron' Rockctdyne (Re£, 6- 13), indicate that con- 
dcni:aiion of tbe exhaust ^ases along tne internal v.;all can occur when 
chc engine is operating .at reduced- power levels. Figure 6-27 shows 
the predicted engine internal wall temperature when the engine is 
operating at the normal, emergency,' and minimum power levels. 

Also shown in this figure are the cor responding saturation tempera- 
tures of the I'R^O er-ihaust gases (i. e. , temperature at which condensa- 
tion of the gas will start), • ■ * 

Since -the v/ali temperature during noimal and emergency power 
G 2 :^oraticn is above the saturation tenrpcr-atur’e, no condensation would 
be anticipated. However, as shown in Figure*. 6-27, the wall tempera- 
ture during minimum power level operation is below the saturation 
temperature and, therefore, condensation of the exhaust, gases in the 
engine wall boundary layer w^ould be expected. The significance of the 
condensation is its tendency to reduce the base region gas recovery . 
temnerature and convective heat transfer. 
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6. 5 


FHOW REVERSAH PHASE ANALYTICAL MODE 


An analytical model of multi-engine flow z^eversal was developed 
and evaluated as a part of this study. Details of the computer program . 
and imodification are discussed in VoTume III of this report. Included 
in Vclunre III are paranretric studies of the Saturn base heat transfer 
and comparisons between scale model test results and the flow reversal 
Q T- o g r a nr . 

Shown in Figures 6-28 through 6-36 is a comparison of the flight 
and scale model measured parameters and the analytical results obtained 
fronr the flow reversal program. 
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Figure 6-23 shows the pressure ir-easurcd in the Saturn flame 
shield region* Also shown in the figure is the analytically computed., 
chching pressure. Chokine pressure is the pressure at which the 


seen in Figure 6-28, the choking pressure measured during the flight 
v/as apprcxi'rnately 2^7 psia. With one of the four center engines cutoff 
the choking pressure was reduced to approximately 0-7 psia. The chok- 
ing pressure co-naputed, using the flow reversal program, is approximately 
0* 24 psia. The analytical choking pressure was computed using a Nash 
factor of unity. The Nash factor controls the anaount of boundary layer 


Lima wnicn is ass- 


id to be reversed. A Nash factor of unity essentially 


: sumes z 


hat the maiArnurn amount of fluid which can be reversed will be 


iversed- Likewise, -a Nash factor of less than unity assumes that 'some 


fluid which might reverse will not be reversed. 

The cQ-mputed choking pressure is significantly less than the 


mcasur 


red choking pressure (Figure 6-28). This appears to be a com- 


mon char acceristic of the flow reversal urogram and, at present, can- 




ed; however, several possibilities can be suggested. Since 


computed pressure is lower than the measured pressure, - this indi- 


car os mar mo: 


fluid is actually reversed during flight than that which 


is being convputed. In Section 6. 2, it was shown that the actual plume 
expandod naore upon leaving the nozzle' than was inidicated by the inviscid 
Low approximations (Figure 6-7). This suggests that either the inviscid 




f the plume boundary is not accurate for large liquid engine 


or a rmexer iriixm': 


layer exists along the outer periphery of the plume. 


Mthcr of these conditions would cause more gas to be reversed into the 


Paso region resulting in a Ligner cnoxing pressure. 
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6-29 shows the ^neasured and computed Saturn S-I flame 
sliiefd vas semocrature. T'ne computed values shov.m are for Nash fac- 
tors ci i.O and 0,35. In Figure 6-29 the computed gas temperature 
acrrees v/cll with the measu^'6d gas teinperature; however/ this is not 
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ic cf the program. As will be seen in later comparisons. 
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tures tend to be hio'her than the imeasured 


:<jrnncrarurcs. 


, The Saturn S-IC flight and scale model measured base gas 
:im ceraiurc s and the computed gas temperatures are shown in Figure 
-30, The Flight temperatures' are those measured by the heat shield 


rct;C; in. 3 men 


fiov/ the surface) C52-106. Because’ of the scale of 


the mode/ different recovery temperatures are obtained from the flow 
reversal program for the model and the flight vehicle. In Figure 6-2^’ 
die computed gas tc-mperature is approximately twice as high as the 
fli‘7ht measured tcrmocraturc. The aas temperature computed for the 
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-ica.lc ‘model shows bettor agreement. As discussed in Section 3, the 
flight gas tom-pcratux*c probes had the potential for error, which naight 


c>:*Diam a pornon oi ‘cno clmc 


between the coimouted temperatures 


ind.'dhc -measured values; Iiovucvcr, it is unlikely that the full difference 


CQxnci DC ■attrioutca to prooe errors. 


* I 1 r/ y S ^ • e 


nodei, and analytical now reversal program convective 


heat flux. The comiP'uted results show good agreement with both the 
flierht and scale ‘model measured heatixi^ rates.* Both the magnitude of 


tne near siux an: 


cho trend predicted for the scale 'fnedei, with respect 


to the fhlght configuration (i. e. , scale model convective heating higher 
•:ha‘n the flight co‘nvective heat flu>:), show good agreement. Apparently, 


the high, predicted gas temperature and the low predicted pressure tend 
to cormoe‘nsate each other when- the convective flux is computed. 


CD CD 


) 
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leisure is t:ie corr.parisor. oi tne :^aturn upper 

:c auaiyrical flow reversal program computed heat shield 
he -orccicted choking pressures indicated choked flow 


altitude -orior to the altitude at which S-II fli 


gnt 3 rage 


L:on occurrec 


olv means that the flow in the S-II 




is predicted to choke immediately after engine ignition and 


u:n:a:n encxea tnrougnour rne ilignr. As can oe seen in rigure o-od, 
compmed base pressure is approximately the same for both the 
acdel and flight configurations. The computed pressures are signi- 
icantly lower than those measured during the flight. During the flight 
ha measured pressure varied betv/een 0. 05 psia and 0. 06 psia. The 
rressure commuted using the flow reversal program varied between 


N 1 and 0, 001 for N = 0. 


• Une cizierence cannor oe 


--rp^aineu at this time except for tne reasons suggested for* tne Saturn* 
flame shield. 

The flight measured and computed S~II' stage base gas tempera- 
izm-s are sriown in Figure 6-33. As with the S-IC stage, the computed 
cmpvjratures are approximately twice as high as the. measured gas 
emneratur es. 


Shown in Figure 6-34 is a comparison of the S-II stage flight, 
model m.easured, and comiputed convective heating rates. As 
ihe S-IC stage the computed heating rates appear to shov/ reasonabh 
: agreement with both the scale mLodel and flight/test results. 

Fi^.ures 6-35 and 6-36 show comparisons, betv/een the S-IY stage 
nlg-iv. measured and the flow reversal model computed heating rates 
pressures. These results indicate approximately the same trend 
blished by the previous conrparison. 
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^ DIAGXCSTIC AND 


Some cc:-clusior.$ arrived at regarding the diagnostic and predic- 
>r, tech.nicues arc the foilov/ing: 


O - ' * ^ V. C' J i 1 - C V 




durin^- the 




as inipox'tarir as riae otner sources ox iioat 
appears that the base convective cooling 
asoiratioix nhasc can be 'ox'cdictcd to an 


n 


: t 


The inability to predict the actual plume boundary 
(inchacling the viscous mixing arid reacting iayex*) 
prevents establishing the onset of flow reversal. A 
thicker and' less energetic plume boundary layer 
co'ald also explshn some of the apparent problems 
witi'j. the analytical 2^-ediction of flow X'eversal. • 


The use of gaseous radiation techniques and band 
modclrj will probably eventually lead to appropriate 
methods for predicting plume -radiation. ' However, 
at the resent tin-xe, these n-xetheds, aiopear to be 
c o mp lex, .t 1 m e - c o n s u nr i n g , and 1 i mii ted ( e x c c p t f o r 
S'pecial comnxataticn) to single loiumes* The latter is 
related *:o tire inability to comToietely die fine the very 
ccmpolicated irnminging plunre thernrodynamic proy^- 
ertres axrd not to the radiation computation teclnriques, 
per se. 

Tor engineez'ing com'puta.iions^y the use of the radiation 
decay method appears to give acceptable results; how- 
evexp this method requires that the plume ground level 
radia'croir oe Known. 

Studies of engine cryopropeilazrt effects indicate that 
ccndcnsarlon of the engine erdraust products can occur 
whexr dre wail tem*oerature is below the saturation 


lire eitec 


oz cenaen- 


ternperazure oi tire exiraust gas. Tire eifects oz cena 
sation on the base environment could be sienificant. 


o The axraiytical flow reversed program appears to over- 
predict the base gas te'crperature and underpredict the 
base pressure; however, these, effects appear to be 
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set v/hc:: predicting the base coii\ 
*. *Tr.e ■program anoears tO oxxer 
making accepiabie pr edictions of 


active heat trails- ' 
an alternate 'means 
the convective heat 


insicr during flow reversal.’ 



of the flew reversal program to accurately 
IS temperature and pressure appears to be 
ufficient reversed mass from the jet* 

, could be related to the inability to accu-’ 


rascly predict the expanded plume mixing layer. 
Further .research in this area is recommended. 
Fmphxisis should be placed on defining the thermo- 
(i/naxnlc and flow fields’ in the mixing uones of jet 
plumes of large rocket engines.. The effe cts of turbine 
exhaust gas discharge and cold walls on the nozzle 
boundary layers should be included'. ■ , 
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